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Abstract 


Micromechanical  crystal  plasticity  finite  element  simulations  of  the  response  of 
synthetic  titanium  microstructures  are  carried  out  with  the  goal  of  quantifying  the 
effect  of  microstructure  on  mechanical  properties.  Two  separate  materials  are  stud¬ 
ied:  (1)  an  a-f3  Ti-6A1-4V  material  and  (2)  a  highly-textured,  rolled  a  Ti-3A1-2.5V 
sheet  material.  Performing  accurate  finite  element  analyses  begins  with  accurate 
image-based  characterization  of  the  morphological  and  crystallographic  features  of 
the  materials  at  the  microstructural  scale.  Then,  statistically  equivalent  represen¬ 
tative  3D  microstructures  are  built  and  meshes  are  generated  for  crystal  plasticity 
based  finite  element  method  (CPFEM)  analysis.  For  the  Ti-3A1-2.5V  material,  ex¬ 
perimental  results  from  the  displacement  controlled  mechanical  testing  of  dog  bone 
shaped,  rolled  specimens  are  used  for  the  calibration  of  elastic  parameters  as  well  as 
anisotropic  crystal  plasticity  parameters.  The  inspection  of  micrographs  of  the  rolled 
material  showed  elongated  grain  shapes  which  led  to  the  updating  of  the  crystal  plas¬ 
ticity  model  to  include  grain  aspect  ratio  dependence  on  the  Hall-Petch  size  effect-an 
update  of  a  previous  size  effect  model  which  assumed  spherical  grains.  Model  vali¬ 
dation  is  achieved  by  comparing  load  controlled  experimental  results  with  simulated 
creep  results.  For  the  Ti-6A1-4V  material,  the  robust  and  validated  analysis  tool  is 
used  to  perform  sensitivity  analyses  and  a  quantitative  understanding  of  how  indi¬ 
vidual  microstructural  parameters  affect  the  mechanical  response  properties  of  the 


alloy  is  developed.  Functional  dependencies  are  proposed  that  directly  connect  the 
metal’s  microstrnctnral  features  to  creep  response,  yield  strength  response,  and  tensile 


response. 
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Chapter  1:  Introduction 


Titanium  alloys  are  widely  used  in  a  number  of  applications  spanning  a  range  of 
industries  including  the  aerospace,  medical,  and  even  sporting  goods  industries  [25]. 
The  positive  properties  of  titanium,  including  high  strength  to  weight  ratio,  good 
fracture  toughness,  and  good  corrosion  resistance  make  this  material  desirable.  One 
downside  to  titanium  alloys,  however,  is  that  they  exhibit  a  ’’cold”  creep  phenomena 
where  time-dependent  deformation,  which  is  normally  observed  at  elevated  temper¬ 
atures  in  most  alloys,  is  seen  to  occur  at  room  temperature  and  at  loads  as  low  as 
60%  of  yield  strength  [38,  39,  28,  20,  29,  35,  43,  15].  The  ’’cold”  creep  mechanism, 
being  at  low  homologous  temperatures,  is  not  diffusion-mediated,  but  occurs  due  to 
dislocation  glide  where  dislocations  are  arrayed  in  a  planar  fashion.  This  planarity  of 
slip  has  been  attributed  to  the  effect  of  short  range  ordering  of  Ti  and  A1  atoms  on 
the  hep  lattice  [29]. 

More  recent  research  has  shown  the  ’’cold”  creep  response  is  strongly  dependent 
on  crystallographic  orientation  in  titanium  [38,  1].  Ghosh  et  al.  have  developed 
CPFEM  models  for  analyzing  the  time-dependent  mechanical  response  of  representa¬ 
tive  titanium  microstructures  [14,  7,  42],  These  models  are  capable  of  capturing  the 
anisotropic  creep  response  of  realistic  titanium  microstructures  with  a  large  number 
of  grains  (on  the  order  of  1000  grains).  In  particular,  the  most  recent  models  are 
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able  to  treat  the  more  complex  transformed  f3  colony  regions  consisting  of  alternating 
layers  of  both  hep  and  bcc  phases.  One  goal  of  this  current  work  is  to  study  the 
anisotropic  creep  response  of  these  relatively  complicated  bimodal  titanium  alloys. 

It  is  well  known  that  the  microstructural  features  of  polycrystalline  materials  di¬ 
rectly  affect  mechanical  properties  (including  the  creep  response  described  above). 
However,  much  of  the  understanding  of  the  relationships  between  microstructure  and 
properties  is  qualitative.  One  approach  to  developing  more  quantitative  understand¬ 
ing  is  through  generating  large  quantities  of  different  microstructures  through  dif¬ 
ferent  heat  treatments  and  thermomechanical  processing.  Then,  through  testing  of 
these  samples,  large  databases  of  experimental  results  are  populated  and,  through 
neural  network  fuzzy  logic  modeling,  functional  dependency  relationships  are  created 
[40,  6,  23].  In  addition  to  the  large  cost  associated  with  processing  and  testing  real 
samples,  one  of  the  downsides  to  this  experimental  approach  is  that  it  is  very  dif¬ 
ficult  to  vary  one  microstructural  feature  while  holding  the  others  constant  during 
heat  treatment  processing  because  the  evolution  of  most  of  these  important  features 
is  highly  coupled. 

Another  method  for  determining  microstructure-property  relationships  is  through 
sensitivity  studies  performed  using  crystal  plasticity  based  finite  element  method 
(CPFEM)  analyses  of  synthetic  microstructures.  This  method  is  pursued  in  the  cur¬ 
rent  work.  One  of  the  beneficial  aspects  of  this  approach  is  that  a  large  array  of 
synthetic  samples  with  different  microstructural  features  can  be  generated  and  tested 
at  a  low  cost  compared  with  testing  of  real  specimen.  Furthermore,  in  contrast  to 
the  limitations  when  generating  real  samples,  it  is  possible  to  vary  the  microstruc¬ 
tural  features  of  these  synthetic  samples  independently  in  order  to  gain  a  more  direct 
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understanding  of  their  effect  on  properties.  This  approach  is  not  without  its  own 
difficulties  and  drawbacks,  which  include  the  need  for  calibration  and  validation  of 
the  models  and  large  amounts  of  computing  power.  However,  there  have  been  recent 
advances  in  the  validation  and  robustness  of  CPFEM  models  (as  discussed  above) 
and  one  goal  of  the  current  work  (see  Chapter  3)  is  to  add  to  the  robustness  of  the 
calibration  and  validation  of  these  models.  Furthermore,  the  increased  availability 
of  the  computing  power  necessary  to  solve  high  resolution  CPFEM  problems  with  a 
large  number  of  grains  means  computational  simulation  has  a  lot  to  offer  towards  the 
solution  to  this  problem.  In  [41],  the  sensitivity  of  yield  strength  to  phase  volume 
fraction  was  studied  using  synthetic  microstructures  consisting  of  cubic  grains.  In 
the  current  work,  more  realistic  grain  morphologies  are  used  and  a  fuller  set  of  re¬ 
sponse  parameters  has  been  studied  including  not  just  yield  strength,  but  also  creep 
parameters  and  Ramb erg- Osgood  parameters. 

The  process  for  developing  a  quantitative  understanding  of  microstructure-property 
relationships  using  CPFEM  analyses  begins  with  the  detailed  characterization  of 
microstructural  features  of  real  titanium  samples.  This  process  consists  of  2D  mi¬ 
crostructure  characterization,  3D  microstructure  model  creation,  meshing  and  CPFEM 
Post-processing  of  the  CPFEM  results  give  stress-strain  curves  and  creep  curves 
whose  parameters  can  be  determined  and  then  studied  with  respect  to  changes  in 
microstructural  features.  The  quantitative  characterization  of  material  microstruc¬ 
ture  is  the  crucial  first  step  in  this  process.  This  task  is  accomplished  through  SEM 
image  analysis,  EBSD  data  collection,  stereological  procedures  and  statistical  quan¬ 
tification.  Once  the  characterization  data  sets  are  acquired,  a  3D  reconstruction  code 
described  in  [8]  is  used  to  create  a  statistically  equivalent,  synthetic  microstructure 


3 


which  is  meshed  and  used  in  computational  simulation.  Results  are  compared  with 
experimental  test  results  and  the  model  is  validated.  Finally,  an  array  of  synthetic 
microstructures  with  varying  statistics  is  generated  and  analyzed  to  gain  an  under¬ 
standing  of  response  sensitivities.  See  Figure  1.1  for  a  flowchart  that  portrays  the  full 
methodology  and  process  described  in  the  current  work. 

Various  aspects  of  this  process  were  carried  out  using  two  separate  but  related 
materials  and  the  work  with  these  two  materials  is  described  in  the  following  2  chap¬ 
ters:  Chapter  2  describes  work  based  on  the  a-/3  Ti-6A1-4V  material,  and  Chapter  3 
describes  work  done  based  on  the  a  Ti-3A1-2.5V  material. 
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Figure  1.1:  Flowchart  of  modeling  methodology  and  process. 
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Chapter  2:  Ti-6A1-4V  Material 


This  chapter  describes  work  done  based  on  the  a-/3  Ti-6A1-4V  material.  Mi¬ 
crograph  images  and  EBSD  scans  of  the  pre-tested  material  and  a  yield  strength 
experimental  data  point  from  a  mechanical  test  of  a  real  sample  provided  a  basis 
for  developing  the  computational  model  and  proposed  functional  forms  connecting 
microstructure  to  mechanical  properties.  The  proposed  model  which  connects  a  lath 
thickness  to  yield  strength  is  validated  against  a  previously  published  linear  regres¬ 
sion  model  generated  from  the  testing  of  75  separate  Ti-6A1-4V  samples  with  various 
heat  treatment  histories. 

2.1  CPFEM  Simulation  of  Polycrystalline  Ti-6A1-4V 

In  order  to  accurately  simulate  the  micromechanical  response  of  polycrystalline 
materials,  it  is  important  to  develop  both  a  good  microstructure  model  as  well  as  a 
good  micromechanical  model.  These  separate  but  related  aspects  to  the  simulation 
are  discussed  in  the  following  pages  beginning  with  the  microstructural  simulation 
model. 
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2.1.1  3D  Polycrystalline  Microstructural  Simulation  and  Mesh 
Generation 

The  following  section  discusses  the  2D  characterization  and  generation  of  3D  mi¬ 
crostructure  models  from  the  2D  data.  It  has  been  shown  that  3D  microstructures 
are  able  to  be  collected  by  a  number  of  techniques,  such  as  Focused  Ion  Beam  (FIB) 
sectioning  [10],  manual  polishing  [33]  and  x-ray  methods  [22],  but  all  remain  some¬ 
what  tedious  and  time  consuming  and  have  not  progressed  to  a  point  where  they  are 
feasible  in  this  type  of  sensitivity  study.  Rather,  the  work  presented  here  will  employ 
an  extrapolation  process  of  estimating  3D  statistics  from  2D  measurements.  This 
area,  known  as  stereology,  has  been  a  source  of  ongoing  investigation  for  years  [34], 
and  will  serve  in  conjunction  with  more  recent  direct  3D  techniques,  to  inform  the 
process  used  in  this  work. 

Microstructural  Characterization  of  the  a-(3  Ti-6A1-4V 

To  complete  the  microstructural  characterization  of  the  a-/3  Ti-6A1-4V,  micro¬ 
graphs  of  the  Ti-6A1-4V  material  were  analyzed  and  statistics  were  compiled.  One 
of  the  micrographs  analyzed  is  shown  in  Figure  2.1.  The  microstructure  consists  of  2 
major  regions:  (1)  transformed  f3  colonies  consisting  of  alternating  a  laths  (hep)  and 
ribs  (bcc)  and  (2)  equiaxed  primary  a  grains  (hep). 

Schematics  of  the  hep  and  bcc  unit  cell  crystal  structure  and  slip  systems  based 
on  atomic  packing  are  shown  in  Figures  2.2  and  2.3.  The  hep  crystal  has  an  non- 
orthogonal  basis  while  the  bcc  crystal  has  an  orthogonal  basis. 

The  crystallographic  orientations  of  this  material  are  influenced  by  the  /3-to-a 
phase  transformation,  which  generally  obeys  the  Burger’s  orientation  relationship  and 
the  resultant  material  has  a  misorientation  distribution  with  preferences  to  certain 


7 


Figure  2.1:  One  of  7  images  at  1000X  magnification. 

variants  [32],  Figure  2.6  shows  the  misorientation  distribution  (MoDF)  of  the  sample 
used  in  this  work.  The  MoDF  is  comparable  to  that  seen  in  the  Ti-6A1-4V  work  done 
by  Randle  in  2008  [32], 

The  sample  was  imaged  using  an  FE1  Sirion  SEM.  A  number  of  images  were 
acquired  using  this  microscope;  including  three  sets  of  seven  images  at  magnifications 
of  750X,  1000X,  and  1500X.  Figure  2.1  shows  one  of  the  images  taken  at  1000X 
magnification.  These  images  were  analyzed  and,  using  the  stereological  techniques 
and  assumptions  described  in  [36]  the  data  shown  in  Table  2.1  were  collected.  These 
particular  morphological  features  are  considered  by  experimentalists  to  be  among  the 
most  important  features  of  these  alloys  that  directly  affect  the  response  [23,  18]. 
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Figure  2.2:  Hep  crystal  structure  and  slip  systems  with  non-orthogonal  basis.  Figure 
taken  from  Balasubramanian  [1], 


Figure  2.3:  Bcc  crystal  structure  and  slip  systems  with  orthogonal  basis.  Figure  taken 
from  Venkatramani  [42], 
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Par  am. 

Description 

Value 

D 

average  grain/col.  size  (ESD) 

11.9  /ini 

Ssd 

stand,  dev.  grain/col.  size 

5.22  fim 

la 

a  lath  thickness 

0.29  /irn 

b 

f3  rib  thickness 

0.089  /.irn 

vf 

vol.  fract.  of  glob,  a  phase 

49.9  % 

VfT 

vol.  fract.  of  total  a  phase 

93  % 

Table  2.1:  Morphological  data  determined  from  SEM  images  and  stereology. 


Grain/colony  size  is  defined  by  the  equivalent  sphere  diameter,  ESD.  ESD  is 
determined  by  first  calculating  equivalent  circle  diameters  ( ECD )  using  image  anal¬ 
ysis  techniques  along  with  the  assumption  that  the  grains/colonies  are  spherical  in 
3D  space  and  therefore  their  size  can  be  approximated  by  the  following  stereological 
formula  which  connects  the  ECD  of  a  sample  with  the  ESD  of  a  sample: 

~ESD  =  -EJCD.  (2.1) 

7T 

According  to  [9],  the  lognormal  standard  probability  density  function  provides  a 
reasonable  fit  to  the  grain/colony  size  data  and  this  was  confirmed  for  the  current 
work  by  the  testing  of  approximately  50  standard  distributions  using  a  maximum 
likelihood  estimation  approach.  This  function  is  defined  by  two  parameters,  average 
and  standard  deviation.  For  these  reasons,  the  lognormal  standard  distribution  will 
be  used  in  this  work  (see  Figure  2.4)  and  will  be  sampled  in  the  synthetic  microstruc¬ 
ture  generation  technique  described  in  the  3D  reconstruction  section  of  this  paper  to 
produce  statistically  equivalent  synthetic  structures. 

In  addition  to  the  above  morphological  analysis,  the  sample  of  Ti-6A1-4V  was  an¬ 
alyzed  crystallographically.  Crystallographic  orientation  data  was  obtained  through 
an  EBSD  scan  containing  approximately  1000  grains/colonies.  The  scan  was  obtained 
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using  an  FEI  Quanta  SEM  and  processed  using  the  code  described  in  [8]  to  provide  hep 
orientation  data  as  well  as  misorientation  and  microtexture  distributions  (see  Figures 
2.5,  2.6,  and  2.7  for  plots  of  the  data).  Misorientation  is  defined  from  the  2D  EBSD 
scan  along  the  grain  boundary  length  but  assumed  as  equivalent  to  a  grain  boundary 
surface  area  in  the  3D  synthetic  structure  generation.  This  assumption  is  reasonable 
in  that  matching  is  done  based  on  a  unitless  fraction  obtained  by  normalizing  the  data 
by  either  total  grain  boundary  length  (in  2D)  or  total  grain  boundary  surface  area 
(in  3D).  The  orientations  of  the  /3  ribs  (which  are  not  directly  obtained  by  the  scan 
because  of  their  small  size)  are  uniquely  defined  from  the  Burger’s  relationship  by  the 
orientation  of  the  a  phase  and  the  lamellar  structure  of  the  transformed  f3  colonies 
and  are  incorporated  into  the  CPFEM  model  based  on  a  mixture  rule  described  in 
[16], 

Microstructural  Simulation  Procedure 

The  microstructural  simulation  is  performed  using  the  code  described  in  [8]  where 
statistically  equivalent  microstructures  are  generated  whose  morphological  and  crys¬ 
tallographic  statistics  are  matched.  Although  the  grains/colonies  have  been  assumed 
to  be  spheres  for  the  purpose  of  characterization,  when  it  comes  to  generation  of 
the  3D  microstructures  the  grains/colonies  are  allowed  to  perturb  slightly  off  sphere 
with  no  correlation  between  the  size  and  shape  of  grains/colonies.  The  morphologi¬ 
cal  orientation  of  the  slightly  elongated  grains/colonies  was  assumed  to  be  random. 
Grains/colonies  are  placed  in  the  synthetic  structure  based  on  neighborhood  con¬ 
straints  assuming  an  average  number  of  neighbors  to  be  approximately  14  grains 
with  a  variation  according  to  grain  size.  This  choice  was  based  on  the  analysis  of 
IN100  in  [9]  in  3D,  where  it  was  seen  that  the  number  of  neighbors  of  a  grain/colony 
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correlated  strongly  to  its  size.  This  assumption  of  correlation  between  the  number  of 
neighbors  of  a  grain/colony  to  its  size  implies  a  lack  of  clustering  of  similarly  sized 
grains/colonies  (i.e.  random  neighborhoods),  which  appears  to  be  valid  when  viewing 
the  2D  micrographs.  After  the  morphological  voxelized  structure  has  been  built,  the 
grains/colonies  are  assigned  hep  orientations  based  on  a  random  sampling  from  the 
ODF.  Misorientation  and  microtexturing  statistics  are  matched  by  an  iterative  pro¬ 
cess  where  orientations  are  allowed  to  switch  between  grains/colonies  or  be  replaced 
by  new  random  orientations  all  the  while  the  error  is  tracked  and  compared  to  sample 
statistics  until  convergence  is  attained. 

Microstructural  Simulation  Validation 


The  microstructural  simulation  procedure  is  validated  through  comparing  the  sam¬ 
ple  statistics  with  the  statistics  of  the  simulated  microstructure.  Additional  valida¬ 
tion  is  achieved  through  a  convergence  study  of  the  algorithm.  See  Figures  2.4,  2.5, 
2.6,  and  2.7  for  graphical  comparisons  of  the  sample  statistics  with  a  500-grain  syn¬ 
thetically  generated  3D  microstructure.  Good  agreement  allows  us  to  consider  the 
synthetic  structure  to  be  considered  statistically  equivalent. 

The  synthetic  microstructure  generation  code  can  be  used  to  generate  structures 
with  any  number  of  grains/colonies.  Figures  2.4,  2.5,  2.6,  and  2.7  show  statistics 
of  a  500-grain/colony  structure,  ft  has  been  seen  that  by  500  grains/colonies  these 
statistics  have  converged  to  a  small  error  as  defined  by  root  mean  square  errors  for 
the  average  and  standard  deviation  of  grain/colony  size.  The  root  mean  square  error 
for  average  grain/colony  size  is  defined  as: 


E 


av 


(2.2) 
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Figure  2.4:  Grain/colony  size,  D,  distribution  comparison.  Nf  is  the  number  fraction 
of  grains/colonies. 


where  A™  is  the  average  grain/colony  size  of  structures  generated  with  %  grains/colonies 
and  ESD  in  this  case  equals  11.9  /mi. 

The  root  mean  square  error  for  the  standard  deviation  of  grain/colony  size  is 
defined  as: 


Esd  — 


(2.3) 


where  Asd  is  the  standard  deviation  of  grain/colony  size  of  structures  generated  with 
i  grains/colonies  and  Ssd  is  the  actual  standard  deviation  calculated  from  the  data 
taken  from  the  EBSD  scan  which  in  this  case  is  5.22  /mi.  Convergence  plots  for 
the  synthetic  structure  generation  algorithm  show  convergence  for  the  matching  of 
the  grain/colony  size  statistics  of  average  and  standard  deviation  (see  Figure  2.8). 
Synthetic  structures  of  less  than  approximately  300  grains/colonies  are  seen  to  not 
match  the  statistics  very  well.  Structures  of  1000  grains/colonies  and  above  show 
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(a)  (b) 


Figure  2.5:  a  phase  (0002)  pole  figures  for  (a)  Sample  data  and  (b)  Synthetic  structure 
as  well  as  (c)  pole  figure  point  density  (PFPD)  distributions  of  the  2  pole  figures  for 
a  more  quantitative  comparison  of  the  crystallographic  matching. 
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Figure  2.6:  Misorientation  distribution  (MoDF)  comparison.  GBf  corresponds  to 
either  grain/colony  boundary  length  fraction  or  grain/colony  boundary  surface  area 
fraction. 


NfN  <15° 


Figure  2.7:  Microtexture  distribution  comparison  where  microtexture  is  defined  by 
the  number  fraction  of  neighbors  (NfN)  with  misorientation  less  than  15°.  Nf  is  the 
number  fraction  of  grains/colonies. 
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only  slight  improvement  and  because  larger  structures  become  computationally  very 
expensive,  structures  of  between  500-600  grains/colonies  were  chosen  for  the  current 
work. 


(a) 

Figure  2.8:  Convergence  Plots:  Error  of 
deviation,  ES(]  (b)  of  grain/colony  size  vs 
microstructure. 


/ 

(b) 

average,  Eav  (a)  and  Error  of  standard 
No.  of  grains/colonies  ( % )  in  simulated 


Mesh  Generation 

The  resultant  structure  of  the  3D  microstructure  model  generation  process  is  a 
voxelized  volume  with  individual  grains/colonies  having  a  phase  identification  and 
an  orientation  defined  by  3  Euler  angles  (see  Figure  2.9).  This  voxelized  volume 
needs  to  be  meshed  for  finite  element  analysis.  While  this  structure  could  be  directly 
transferred  into  the  finite  element  environment,  the  voxelized  nature  of  the  structure 
causes  a  ’’stair-stepped”  boundary  between  grains/colonies.  The  ’’jagged”  boundaries 
have  been  shown  to  be  the  source  of  local  artifacts  during  simulations  [17].  The 
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current  work  uses  the  4-noded  TET4  element  instead  in  an  attempt  to  remove  the 
artificial  voxelized  boundaries. 

As  this  element  uses  linear  interpolation  functions,  the  strains  are  constant  through¬ 
out  the  element  and  the  element  is  integrated  numerically  using  one  point  at  the 
centroid  of  the  element.  Even  though  the  shape  functions  for  displacement /strain 
are  linear,  for  the  500  to  600-grain/colony  microstructures  used  for  the  current  work 
the  authors  have  selected  meshes  that  are  refined  enough  to  the  end  that  the  creep 
and  constant  strain  rate  responses  have  converged.  In  other  words,  there  is  little 
sensitivity  of  further  mesh  refinement  to  response. 

The  commercial  meshing  software,  Simmetrix  [37],  is  used  to  generate  the  finite 
element  mesh  (see  Figure  2.9).  First,  a  triangular  surface  mesh  is  generated  along  the 
interior  grain  boundaries  and  cube  boundaries.  Then,  this  triangular  surface  element 
mesh  is  extended  into  the  full  3D  volumetric  tetrahedral  mesh.  The  meshes  used 
for  the  current  work  contain  between  approximately  100,000  to  120,000  elements  and 
20,000  to  22,000  nodes. 

The  meshes  are  checked  for  distorted  elements  and  it  has  been  seen  that  the 
meshes  used  for  the  current  work  have  a  very  small  number  of  elements  with  aspect 
ratio  of  40  or  higher  (on  the  order  of  0.01%).  Because  the  current  work  is  concerned 
with  the  volume  average  creep  and  constant  strain  rate  responses,  it  is  reasonable  to 
concede  this  very  small  amount  of  distorted  elements. 

Two  different  microstructures  were  generated  that  are  statistically  equivalent:  one 
with  500  grains/colonies  and  one  with  600  grains/colonies.  Increasing  the  number  of 
grains/colonies  in  the  simulated  microstructures  above  500  shows  no  affect  upon  the 
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Figure  2.9:  Voxelized  volume  (a)  before  and  (b)  after  meshing  for  the  600-grain 
microstructure.  The  cube  length  dimension,  l0 ,  is  68  /jm.  The  grayscale  contour 
variable,  C  or.,  is  the  c-axis  orientation  defined  as  the  angle  between  the  c-axis  of  the 
hep  phase  of  the  grain/colony  and  the  loading  direction  which  in  this  case  is  2. 


simulated  micromechanical  creep  response  parameters  A  and  m  (discussed  in  Section 
2.2).  The  model  has  converged  with  respect  to  these  parameters. 

2.1.2  Micromechanical  Analysis  of  Simulated  Polycrystalline 
Microstructures  Using  a  Crystal  Plasticity  Constitu¬ 
tive  Model 

An  isothermal,  size-dependent  and  rate-dependent  crystal  plasticity  finite  element 
computational  model  described  and  developed  in  [14,  7,  42]  is  used  in  conjunction 
with  MSC/Marc  Mentat  and  an  in-house  parallelized  code  to  simulate  the  response 
of  the  synthetically  generated  microstructures  of  the  bimodal  a  —  (3  Ti  6A1-4V  under 
various  loading  and  boundary  conditions. 
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Constitutive  Model 


Crystal  deformation  is  modeled  by  a  combination  of  elastic  stretching  and  rotation 
of  the  crystal  and  plastic  slip  on  the  different  slip  systems.  The  stress-strain  relation 
is  given  in  terms  of  the  second  Piola-Kirchoff  stress  (S  =  detF<Fe^1aFe^T)  and  the 
work  conjugate  Lagrange  Green  strain  tensor  (Ee  =  (1/2)  |Fe2Fe  —  I})  as, 

S  =  C  :  Ee,  (2.4) 

where  C  is  the  fourth-order  anisotropic  elasticity  tensor,  a  is  the  Cauchy  stress  tensor, 
and  Fe  is  the  elastic  deformation  gradient  defined  by  the  relation, 

Fe  =  FFP_1,  detFe  >  0.  (2.5) 

F  represents  the  deformation  gradient  and  Fp  its  plastic  component.  The  incom¬ 
pressibility  constraint  is  given  by  detFp  =  1.  The  flow  rule  describing  the  plastic 
deformation  is  cast  in  terms  of  the  plastic  velocity  gradient, 

a  V™ 

Lp  =  FpFp_1  =  T  =  ^  —  stgn(Ta ),  (2.6) 

a  9 

where  7“  is  the  plastic  shearing  rate,  r“  is  the  resolved  shear  stress,  ga  is  the  slip 
system  deformation  resistance  (all  variables  are  specified  on  the  ceth  slip  system  of 
the  phase),  m  is  a  material  rate  sensitivity  parameter  and  the  Schmid  tensor,  Sq,  is 
expressed  as 

Sq  =  m“  (g)  n“.  (2.7) 
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The  slip  system  deformation  resistance,  ga,  evolves  along  with  the  hardening  rates 
as  described  in  [21,  12]  as 

nslip 

ga  =  hap  \¥\  =  ^2  <fPhP  1^ |  ,  (2.8) 

13=1  /3 

where  hal3  is  the  strain  hardening  rate  due  to  both  self  and  latent  hardening,  h 13  is 
the  self-hardening  rate,  and  gQ/3  is  a  matrix  describing  the  latent  hardening.  For  the 
hep  a  phase,  it  is  assumed  that  the  evolution  of  the  self-hardening  rate  is  governed 
as 

hp  =  h $  1  -  ^  sign  ^1  ,g%  =  g  ,  (2.9) 

where  /if  is  the  initial  hardening  rate,  g@  is  the  saturation  slip  deformation  resistance, 
and  r,  g,  and  n  are  slip  system  hardening  parameters.  A  different  relation  is  used  for 

the  evolution  for  the  bcc  f3  phase  found  in  transformed  f3  colonies  as  follows: 

~fhP-hP\  ^  rtnslip 

hP  =  hPa  +  seeh2  (  —0 - 1  )  7“  | ip\dt,  (2.10) 

Lvr^-ro/  J  p=i 

where  and  are  the  initial  and  asymptotic  hardening  rates,  rf  represent  the 
saturation  value  of  the  shear  stress  when  /if  =  0,  and  7 a  is  a  measure  of  total  plastic 
shear. 


In  continuum  plasticity,  the  dependence  of  the  yield  stress  on  the  grain  size  has 
been  expressed  by  Hall  and  Petch  in  the  1950s  [11,  30].  A  Hall-Petch-type  equation  is 
used  in  the  crystal  plasticity  formulation.  This  equation  relates  the  initial  slip  system 


deformation  resistance,  ga ,  to  a  characteristic  size  as 


9a  =  9o  + 


(2.11) 


where  g^  and  Ka  are  size-effect-related  slip  system  constants  that  refer  to  the  initial 


slip  system  deformation  resistance  and  slope,  respectively,  and  Da  is  the  character¬ 
istic  length  scale  governing  the  size  effect  (see  Section  3.2.2  for  a  grain  aspect  ratio 
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dependence  model  update  on  this  size  effect).  According  to  work  done  in  [42],  the 
transformed  f3  regions  have  D°  values  which  correspond  either  to  the  colony  size, 
the  a  lath  thickness,  or  f3  rib  thickness  depending  on  ease  of  slip  conditions  at  the 
hcp-bcc  interfaces  of  the  colony  lath  structure.  For  primary  a  grains,  the  D°  values 
are  assumed  to  be  the  ESD  of  the  grain. 

Material  Properties 

Elasticity  and  crystal  plasticity  parameters  are  implemented  and  calibrated  from 
a  multi- variable  optimization  method  developed  and  carried  out  in  [13]  and  [7]  where 
a  Ti-6A1  single  crystals  and  a  —  (3  Ti-6242  single  colonies  were  experimentally  tested. 
Tables  2.3  and  2.2  give  the  elastic  parameters  for  both  the  hep  and  bcc  phases. 
Tables  2.4,  2.5,  and  2.6  give  the  calibrated  values  of  the  plastic  parameters  for  both 
the  hep  and  bcc  phases  in  the  transformed  /3  colonies  as  well  as  the  hep  parameters 
for  the  primary  a  grains.  An  implicit  backward  Euler  time-integration  scheme  is  used 
for  the  solution  to  the  dynamic  problem  using  the  commercial  finite  element  code, 
MSC/Marc  Mental  [27]  along  with  the  user-defined  material  routine  UMAT. 


Cij  parameter 

Value  (GPa) 

A'll=C,22 

170.0 

C33 

204.0 

C'i2=C<2l 

98.0 

Cl3 =C<31  =  C23  =  C32 

86.0 

C44 

Cn-C12 

C,55=C<66 

102.0 

Other 

0 

Table  2.2:  Calibrated  elastic  parameters  of  the  transversely  isotropic  tensor  for  the 
hep  phase. 
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Cij  parameter 

Value  (GPa) 

Cll  =  C,22=C,33 

250.21 

/^i  S~1 

^13— ^31— ^23— (-/32— ^12— ^21 

19.0 

C,44=C,55=C,66 

230.65 

Other  Cij 

0 

Table  2.3:  Calibrated  elastic  parameters  of  the  cubic  symmetric  tensor  for  bcc  phase. 


Parameters  for 
Slip  System 

m 

9o 

(MPa) 

4(S'1) 

ho 

hs 

To 

(MPa) 

rs 

(MPa) 

{101} 

0.02 

450.00 

0.0023 

1500.0 

25.0 

500.0 

200.0 

{112}  soft 

0.02 

429.82 

0.0023 

1371.6 

25.0 

500.0 

200.0 

{112}  hard 

0.02 

409.63 

0.0023 

1979.8 

25.0 

500.0 

200.0 

{123}  soft 

0.02 

451.28 

0.0023 

2312.0 

25.0 

500.0 

200.0 

{123}  hard 

0.02 

400.67 

0.0023 

1634.8 

25.0 

500.0 

200.0 

Table  2.4:  Calibrated  parameters  for  the  bcc  slip  systems  in  the  homogenized  trans¬ 
formed  j3  colonies. 


Parameters  for 
Slip  System 

m 

9o 

(MPa) 

ho 

r 

7 

(MPa) 

n 

(MPa) 

Basal  ai 

0.02 

284.0 

0.0023 

1500.0 

0.30 

450.0 

0.14 

Basal 

0.02 

315.0 

0.0023 

2300.0 

0.30 

634.0 

0.10 

Basal  a3 

0.02 

243.0 

0.0023 

8000.0 

0.40 

371.0 

0.05 

Prismatic  ai 

0.02 

240.0 

0.0023 

3450.0 

0.29 

504.0 

0.15 

Prismatic  &2 

0.02 

210.0 

0.0023 

6500.0 

0.20 

583.0 

0.25 

Prismatic  a3 

0.02 

240.0 

0.0023 

3600.0 

0.30 

504.0 

0.15 

Pyr.  <a> 

0.02 

395.0 

0.0023 

100.0 

0.10 

550.0 

0.01 

Pyr.  <c+a> 

0.02 

623.3 

0.0023 

100.0 

0.10 

1650.0 

0.01 

Table  2.5:  Calibrated  parameters  for  the  hep  basal  <a>  and  prismatic  <a>  slip 
systems  in  the  homogenized  transformed  f3  colonies. 
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Parameters  for 
Slip  System 

m 

9o 

(MPa) 

h0 

r 

7 

(MPa) 

n 

(MPa) 

Basal  <a> 

0.02 

284.00 

0.0023 

100.0 

0.1 

450.0 

0.01 

Prism.  <a> 

0.02 

282.24 

0.0023 

100.0 

0.1 

550.0 

0.01 

Pyr.  <a> 

0.02 

395.00 

0.0023 

100.0 

0.1 

550.0 

0.01 

Pyr.  <c+a> 

0.02 

623.30 

0.0023 

100.0 

0.1 

1650.0 

0.01 

Table  2.6:  Calibrated  parameters  for  the  slip  systems  in  the  primary  a  grains. 


Boundary  and  Loading  Conditions 

For  both  the  constant  strain  rate  and  creep  simulations  to  suppress  rigid  body 
modes,  symmetric  constraint  conditions  are  applied  where  the  back  faces  of  the  cube 
are  constrained  as  ui=0  on  the  1-face,  u2=0  on  the  2-face,  and  u3=0  on  the  3-face. 
The  cube  is  free  to  move  in  all  other  directions.  For  constant  strain  rate  simulations, 
one  of  the  outer  faces  is  given  a  constant  strain  rate  displacement  boundary  condition 
of  Ua(t)  =  l0(exp(e^t)  —  1),  where  lQ  is  the  initial  dimension  of  the  cube  and  ecH  is  the 
applied  constant  strain  rate  with  i  corresponding  to  which  direction  the  cube  is  being 
loaded  in.  For  creep  simulations,  a  constant  load  is  applied  to  one  of  the  outer  faces. 

2.2  Comparison  of  the  CPFEM  Model  with  Experimental 
Results 

The  CPFEM  model  results  are  compared  with  the  experimental  data  available  for 
the  tested  specimen.  One  data  point  was  available;  that  being  the  engineering  yield 
strength  (Y 5=821  MPa)  as  defined  by  a  0.2%  strain  offset  method  for  a  tensile  exper¬ 
iment  conducted  at  a  strain  rate  of  10-4s-1.  Post-processing  of  the  solution  gives  the 
engineering  stress-strain  curve  and  the  simulated  yield  strength  (YS)  is  determined 
to  be  833  MPA  (based  on  the  2%  strain  offset  method)  with  a  small  gfi  parameter 
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calibration  (10%  reduction  required)  deemed  reasonable  because  parameters  were  not 
originally  calibrated  for  this  particular  alloy.  This  YS  value  is  compared  with  the 
experimental  value  of  821  MPA.  The  error  is  1.5%  after  the  slight  calibration.  The 
entire  stress-strain  curve  determined  by  the  simulation  can  be  seen  in  Figure  2.10. 
This  result  is  encouraging  considering  the  limited  2D  data  that  was  available  for 
characterization  of  the  material. 


Figure  2.10:  Engineering  stress-strain  response  with  comparison  of  simulated  YS  and 
experimental  Y S.  The  2-direction  is  the  loading  direction. 


Figure  2.11  shows  1  of  the  many  internal  grains  before  and  after  creep  loading  at 
700  MPa  (approx.  83%  of  yield  strength)  for  10,000  sec.  The  contour  variable,  ep,  is 
plastic  strain.  The  deformation  is  scaled  by  a  factor  of  50. 
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Figure  2.11:  A  single  internal  grain  in  its  (a)  undeformed  and  (b)  deformed  state 
after  creep  loading  in  the  2-direction  at  700  MPa  for  10,000  sec.  The  deformation  is 
scaled  by  a  factor  of  50  and  the  contour  variable,  ep,  is  plastic  strain. 


2.3  Microstructure-Dependent  Macroscopic  Models  for  Yield 
Strength,  Creep  and  Constant  Strain  Rate  Behavior 

Much  advance  has  been  made  towards  developing  microstructure-property  rela¬ 
tionships  for  titanium  alloys  through  the  use  of  neural  network  models  [40,  6,  18,  19, 
44] .  However,  the  time  and  cost  that  is  associated  with  building  those  neural  network 
models  can  be  quite  significant.  They  are  expensive  to  develop  in  that  they  require 
the  population  of  a  huge  database  of  experimental  data  to  be  useful.  Furthermore,  in 
order  to  train  the  neural  network,  the  experiments  need  to  be  performed  on  samples 
that  have  microstructural  features  that  span  a  wide  range  of  values.  To  achieve  this 
variation  in  microstructure,  it  is  necessary  to  perform  expensive  highly  controlled 
heat  treatments  on  the  material  before  testing.  Additionally,  because  of  the  strong 
interrelationship  between  microstructural  features  such  as  grain/colony  size  and  a 
lath  size,  it  is  essentially  impossible  through  experimental  means  to  vary  only  one 
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microstructural  feature  in  order  to  investigate  its  effect  on  mechanical  properties  [23]. 
The  computational  tools  and  methods  discussed  and  implemented  here  provide  an 
alternative  which  helps  to  overcome  these  limitations.  Sensitivity  studies  have  been 
performed  in  the  current  work  to  lay  the  groundwork  for  proposed  functional  forms 
that  relate  microstructural  features  to  mechanical  properties. 

2.3.1  Definition  of  Model  Response  Parameters 

Although  the  CPFEM  micromechanical  model  defines  the  deformation  gradient  as 
F  =  FeFp,  for  describing  the  stress-strain  curves,  true  strain  is  additively  decomposed 
into  its  elastic  and  plastic  portions  as  e  =  ee+ep.  It  is  understood  that  this  assumption 
does  not  hold  well  for  values  of  large  strain;  however,  this  inaccuracy  is  deemed 
acceptable  for  the  current  work  where  parameters  are  fit  using  simulation  data  that 
falls  mainly  in  the  smaller  strain  regions.  The  Ramberg-Osgood  equation  [31]  is  given 
as: 

e=|  +  A'(^r,  (2.12) 

where  the  first  term  represents  the  elastic  portion  of  the  strain,  ee  and  the  second 
term  represents  the  plastic  portion  of  the  strain,  ep.  K  and  n  are  two  dimensionless 
material  parameters.  Post  processing  in  crystal  plasticity  distinguishes  the  plastic 
strain  from  the  total  strain  and  curves  are  fit  to  determine  the  K  and  n  material 
response  parameters.  These  parameters  provide  a  good  description  of  the  constant 
strain  rate  response  of  the  material. 

For  describing  the  creep  response,  a  power  creep  law  given  by  Lubahn  and  Felgar 
is  used  [24]  where  again  the  true  strain  is  decomposed  into  its  clastic  and  plastic 
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portions  and  the  plastic  portion  of  the  response  is  given  by 


ep  =  Atm,  (2.13) 

where  t  is  time,  and  A  (with  units  of  1/sec)  and  m  (which  is  dimensionless)  are  mate¬ 
rial  constants  which  are  fit  to  the  creep  simulation  results.  Again,  it  is  acknowledged 
that  the  additive  decomposition  does  not  hold  well  for  large  strain  problems  but  this 
inaccuracy  is  accepted  based  on  the  fact  that  most  of  the  simulation  fitting  is  in  a 
smaller  strain  regime. 

2.3.2  Sensitivity  Studies 

In  this  study,  an  array  of  synthetic  specimens  is  produced  and  studied  and  sen¬ 
sitivities  of  response  to  changes  in  microstructural  features  are  reported.  A  realistic 
reference  microstructure  is  used  and,  one  by  one,  each  feature  is  varied  while  all  other 
features  are  held  constant.  The  current  work  studies  the  sensitivity  of  response  to 
variations  in  a  lath  size,  la,  /3  rib  size,  1/3,  volume  fraction  of  primary  a  grains,  Vf, 
and  average  grain/colony  size,  D.  As  discussed  above,  this  is  a  unique  capability  that 
is  not  available  through  experimental  means.  The  response  values  that  are  studied 
include  yield  strength,  the  Ramberg- Osgood  parameters  K  and  n  (see  Eqn.  2.12), 
and  the  power  creep  law  parameters  A  and  m  (see  Eqn.  2.13).  The  following  figures 
show  sensitivities  of  these  microstructural  features  to  Y'S  (Figures  2.12-2.15),  as  well 
as  power  creep  law  parameters,  A  (Figures  2.16-2.19)  and  m  (Figures  2.20-2.23),  and 
Ramberg-Osgood  parameters,  K  (Figures  2.24-2.27)  and  n  (Figures  2.28-2.31). 
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Figure  2.12:  Sensitivity  of  yield  strength  (YS)  to  changes  in  a  lath  thickness  ( la ). 


Figure  2.13:  Sensitivity  of  yield  strength  (Y S)  to  changes  in  (3  rib  thickness  (Ip). 


The  studied  range  of  a  lath  thicknesses  and  the  corresponding  (5  rib  thicknesses 
from  approximately  0.2  to  1.0  /im  and  0.1  to  0.3  /im  are  realistic  ranges.  Experimen¬ 
tally,  ranges  that  cover  sizes  differing  by  500%  are  easily  attained  in  real  samples  by 
varying  processing  cooling  rates  [23]. 


Figure  2.14:  Sensitivity  of  yield  strength  (Y S)  to  changes  in  volume  fraction  of  pri¬ 
mary  a  ( Vf ). 

Simulation  results  show  that  in  general  changes  to  la  and  Ip  have  a  smaller  effect 
on  response  compared  to  the  larger  effects  seen  from  changes  in  Vf  and  D.  The 
relatively  smaller  affect  of  these  feature  sizes  on  response  can  be  explained  by  the 
fact  that  Vf  and  D  are  relatively  larger  features  and  therefore  it  is  expected  that 
variations  in  these  features  would  have  a  relatively  greater  affect  on  the  response. 

As  seen  in  Figures  2.12  and  2.13  (as  well  as  the  other  figures  which  investigate 
the  sensitivities  of  la  and  Ip),  the  results  for  a  0%  Vf  sample  were  also  investigated. 
The  0%  Vf  cases  correspond  to  a  100%  transformed  (3  colony  sample.  As  would  be 
expected,  there  is  a  greater  sensitivity  to  changes  in  la  and  Ip  for  the  0%  Vf  cases 
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Figure  2.15:  Sensitivity  of  yield  strength  (Y S)  to  changes  in  average  grain/colony 
size  (D). 


compared  to  the  50%  Vf  cases.  This  is  explained  by  the  fact  that  the  primary  a 
grains  have  no  lath  or  rib  structure.  In  other  words,  the  fewer  of  these  pure-ct  grains 
in  the  structure  and  the  more  of  the  colony  structure  the  higher  the  effect  of  changes 
to  the  lath  and  rib  thicknesses. 
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Figure  2.16:  Sensitivity  of  K  parameter  to  changes  in  a  lath  thickness  (la). 
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Figure  2.17:  Sensitivity  of  K  parameter  to  changes  in  (5  rib  thickness  (Ip). 


v, 


Figure  2.18:  Sensitivity  of  K  parameter  to  changes  in  volume  fraction  primary  a  (Vf). 


Figure  2.19:  Sensitivity  of  K  parameter  to  changes  in  grain/colony  size  (D). 
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Figure  2.20:  Sensitivity  of  n  parameter  to  changes  in  a  lath  thickness  (/„). 
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Figure  2.21:  Sensitivity  of  n  parameter  to  changes  in  f3  rib  thickness  {Ip). 


Figure  2.22:  Sensitivity  of  n  parameter  to  changes  in  volume  fraction  primary  a  (V/). 
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Figure  2.23:  Sensitivity  of  n  parameter  to  changes  in  grain/colony  size  (D). 


Focusing  on  the  sensitivity  of  YS  to  the  various  microstructural  features  (Figures 
2.12-2.15),  the  trend  is  consistent  for  all  features  that  as  more  grain/colony  boundary 
or  more  lath/rib  boundary  is  introduced  into  the  structure  the  YS  increases.  This 
trend  is  in  agreement  with  the  understanding  that  boundaries  impede  dislocation 
motion  and  therefore  hinder  the  onset  of  plastic  deformation  and  ultimately  increase 
the  strength  of  a  material. 


Figure  2.24:  Sensitivity  of  A  parameter  to  changes  in  a  lath  thickness  ( la ). 
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Figure  2.25:  Sensitivity  of  A  parameter  to  changes  in  f3  rib  thickness  (Ip). 
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Figure  2.26:  Sensitivity  of  A  parameter  to  changes  in  volume  fraction  primary  a  (Vf). 

It  is  also  worth  taking  note  of  the  sensitivity  of  YS  to  average  grain/colony 
size  (Figure  2.15).  This  non-linear  response  shows  a  Hall-Petch-type  relationship 
where  initially  (with  smaller  grain  sizes)  there  is  a  very  large  sensitivity  to  changes 
in  grain/colony  size,  D,  and  then  as  D  become  larger  and  larger  there  is  a  saturation 
and  there  is  not  as  much  sensitivity  to  increases  in  D  after  it  reaches  around  20  /mi. 
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Figure  2.27:  Sensitivity  of  A  parameter  to  changes  in  grain/colony  size  (D). 


Figure  2.28:  Sensitivity  of  m  parameter  to  changes  in  a  lath  thickness  (la). 


Figure  2.29:  Sensitivity  of  m  parameter  to  changes  in  f3  rib  thickness  (Ip). 


35 


Figure  2.30:  Sensitivity  of  m  parameter  to  changes  in  volume  fraction  primary  a  (Vf). 
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Figure  2.31:  Sensitivity  of  m  parameter  to  changes  in  grain/colony  size  (D). 
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These  results  are  comparable  with  the  experimental  results  seen  by  Yoshimura  and 
Nakahigashi  [45].  They  have  obtained  grain  sizes  as  small  as  1  //m  in  a-/3  Ti  using 
hydrogen  treatments  and  saw  that  YS  results  followed  the  Hall-Petch  relationship 
well  at  these  values.  A  common  two-parameter  version  of  the  Hall-Petch  relationship 
is  given  as, 

YS  =  a0+YL  (2.14) 

where  cto  and  K*  (not  to  be  confused  with  the  K  parameter  of  the  Ramberg-Osgood 
equation  used  in  this  work)  are  macroscopic  size-effect  constants. 

Plotting  the  YS  as  a  function  of  1/D0-5  (Figure  2.32),  a  value  of  <t0  of  760  MPa  and 
K*  of  240  MPa^/jUm  are  obtained  from  a  straight  line  fit.  These  values  are  compared 
to  those  obtained  by  [42]  for  a  synthetic  specimen  of  Ti-6242  of  750  MPa  and  250 
MPa^/yum  respectively. 


1/D05  (pm-05) 


Figure  2.32:  Hall-Petch  two  parameter  model. 


Focusing  our  attention  to  Figures  2.15,  2.19,  2.23,  2.27,  and  2.31,  the  effect  of 
average  grain  size,  H,  on  all  the  response  parameters  can  be  compared.  Notice  that, 
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just  like  for  YS,  the  K  and  n  parameters,  as  well  as  the  creep  parameters,  A  and 
m  show  the  similar  characteristic  of  high  sensitivity  at  low  grain  sizes  and  a  lessened 
sensitivity  as  grain  size  increases.  In  future  work,  it  could  be  examined  if  for  a  very 
large  average  grain  size  the  effect  of  changes  in  alpha  lath  and  beta  rib  thickness 
would  be  heightened. 

2.3.3  Functional  Forms 

Full  functional  forms  that  connect  microstructural  features  and  material  response 
parameter  are  developed  using  the  sensitivity  information  discussed  in  the  previous 
section.  The  forms  are  developed  assuming  that  variational  effects  are  independent 
of  each  other. 

The  yield  strength,  YS,  as  a  function  of  a  lath  thickness,  la ,  (3  rib  thickness,  Ip, 
volume  fraction  of  primary  a,  Vf,  and  average  grain/colony  size,  D,  is  determined  to 
be 

YS(la,  Ip,  Vf,  D)  =  kYS  *  YS(la )  *  YS(lp)  *  YS{Vf )  *  YS(D),  (2.15) 

where  YS(la)  =  -27.4*/Q  +  852,  YS{lp)  =  -147*^  +  857,  YS(Vf)  =  -163*F/  +  930, 
YS(D )  =  968  *  D-0  055,  and  kys  is  a  calibrated  constant  which  equals  1.62E-9. 

A  validation  of  the  proposed  functional  form  connecting  la  to  YS  is  obtained  by 
comparison  with  previously  published  work.  In  2002,  Tiley  [40]  developed  a  linear 
regression  model  based  off  neural  network  modeling  of  experimental  yield  strength 
results  from  approximately  75  different  Ti-6A1-4V  samples  that  had  been  heat  treated 
to  produce  different  a  lath  thicknesses.  Tiley’s  linear  regression  model  obtained  from 
a  large  amount  of  experimental  data  provides  a  good  validation  for  the  proposed 
functional  form  connecting  YS  to  la.  The  linear  regression  model  obtained  from 
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experiments  given  by  Tiley  as,  Y S(la)  =  —15.4  *  la  +  826,  (MPa),  is  compared  with 
the  model  proposed  here  from  simulations  of  Y'S(la)  =  —27.4  *  la  +  852,  (MPa). 
Figure  2.33  gives  a  plot  comparing  the  two  linear  regression  models. 


Ia  (Mm) 


Figure  2.33:  Comparison  of  the  linear  regression  model  obtained  from  experiment 
[40]  and  that  obtained  from  simulation  for  Y S  vs.  la. 


The  Ramb erg- Osgood  parameter,  K,  as  a  function  of  a  lath  thickness,  la,  (3  rib 
thickness,  Ip,  volume  fraction  of  primary  a,  Vf,  and  average  grain/colony  size,  D ,  is 
determined  to  be 

K(la,  Ip,  Vf,  D)  =  kK  *  K(la)  *  K(lp)  *  K(Vf )  *  K(D),  (2.16) 

where  K(la)  =  9.44£20  *  Z"1-97,  K{lp)  =  1.21E2A  *  Z^88,  K(Vf)  =  1.03E13  *  e39'7*^ , 
K(D)  =  3. 897724  *  D~2A8,  and  kx  is  a  calibrated  constant  which  equals  9.22E-67. 

The  Ramberg-Osgood  parameter,  n,  as  a  function  of  a  lath  thickness,  Za,  f3  rib 
thickness,  Ip,  volume  fraction  of  primary  a,  Vf,  and  average  grain/colony  size,  D,  is 
determined  to  be 

n(la,  Ip,  Vf,  D )  =  kn  *  n(la )  *  n(lp )  *  n(Vf)  *  n(D),  (2.17) 
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where  n(la)  =  —0.794  *  lQ  +  11.6,  n(lp)  =  1.82  *  Ip  +  11.3,  n(Vf)  =  7.65  *  e°-no*vf ; 
n(D )  =  13.0  *  f)-0-00536,  and  kn  is  a  calibrated  constant  which  equals  6.24E-4. 

The  creep  parameter,  A,  as  a  function  of  a  lath  thickness,  la,  [3  rib  thickness,  Ip, 
volume  fraction  of  primary  a ,  Vf,  and  average  grain/colony  size,  D,  is  determined  to 
be 

A(la,  Ip,  Vf,  D)  =  kA*  A(la )  *  A{lp)  *  A(Vf)  *  A(D),  (2.18) 

where  A(la)  =  5.95 E  -05  *la  +  1.80 E  -  04,  A(lp)  =  2.66 E  -  04  *  Ip  +  1.71  E  -  04, 
A(Vf)  =  3.50 E  -05  *Vf  +  2.07 E  -  04,  A(D)  =  1.29 E  -  04  *  ln(D)  -  3.57 E  -  05,  and 
kA  is  a  calibrated  constant  which  equals  9.59E10. 

The  creep  parameter,  m,  as  a  function  of  a  lath  thickness,  la,  [3  rib  thickness,  Ip, 
volume  fraction  of  primary  a,  Vf,  and  average  grain/colony  size,  D,  is  determined  to 
be 

m(la,  Ip,  Vf,  D)  =  km*  m(la)  *  m(lp)  *  m(Vf)  *  m(D),  (2.19) 

where  m(la)  =  1.49E  -03  *  la  +  2.04 E  -  01,  m(lp)  =  1A1E  -02  *lp  +  2.04E  -  01, 
m(Vf)  =  9.98 E  —  02  *  Vf  +  1.51  E  —  01,  m(D)  =  2. 4917  —  01  *  D~9- slE-°2,  and  km  is  a 
calibrated  constant  which  equals  1.18E02. 

2.3.4  An  Extended  Investigation  of  Material  Anisotropy 

With  a  validated  and  calibrated  CPFEM  model,  the  material  anisotropy  is  in¬ 
vestigated  by  loading  a  synthetic  microstructure  with  similar  statistics  in  different 
directions.  In  addition  to  the  experimental  loading  direction  (i=2) ,  the  response  of 
the  synthetic  structure  to  loadings  in  the  other  2  orthogonal  directions  (i=l,3)  is 
also  considered.  The  constant  strain  rate  results  shown  in  Figure  2.34  confirm  the 
anisotropy.  This  anisotropy  is  a  direct  result  of  the  non-random  grain  orientations 
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and  the  anisotropy  in  the  parameters  which  govern  the  clastic-plastic  response  of  the 
individual  slip  systems. 


Figure  2.34:  Simulated  engineering  stress-strain  response  for  loadings  in  3  directions 
(i— 1)2,3). 


The  anisotropy  of  the  microstructure  under  different  loading  directions  can  be 
characterized  by  the  softness  index.  The  softness  index  is  a  dimensionless  parameter 
defined  as, 

SI  =  SFS  *  YS/gs0,  (2.20) 

where  SFS  is  the  schmid  factor  of  the  softest  slip  system  of  the  grain/colony  (which 
depends  on  the  loading  direction  and  the  slip  plane  normal  direction  of  the  softest  slip 
system),  YS  is  the  macroscopic  yield  strength  of  the  entire  sample,  and  is  the  size- 
dependent  initial  slip  system  resistance  of  the  softest  slip  system  of  the  grain/colony. 
Figure  2.35  shows  the  distribution  of  softness  index  for  the  three  loading  direction 
cases  (i=l,2,3).  Softness  index  values  of  1  or  higher  indicate  that  the  grain/colony 
will  plastically  deform  when  the  applied  load  reaches  YS  while  values  below  1  indicate 


41 


that  the  grain/colony  will  not  see  any  plastic  deformation  at  YS.  Notice  that  loading 
directions  i=2  and  i=3  have  more  ’’softer”  grains/colonies  compared  to  case  i=l.  This 
explains  the  ” harder”  response  of  loading  direction  i=l  seen  in  Figure  2.34. 


Figure  2.35:  Distribution  of  softness  index  (SI)  for  the  three  loading  directions 
(i=l ,2,3) .  Nf  is  the  number  fraction  of  grains/colonies.  Notice  that  loading  dir. 
i=l  has  fewer  ’’softer”  grains/colonies. 


To  further  investigate  anisotropy,  modified  versions  of  the  Ramberg-Osgood  equa¬ 
tion  2.12  and  power  creep  law  equation  2.13  are  used  to  characterized  the  anisotropy 
in  the  constant  strain  rate  and  creep  response.  The  Ramberg-Osgood  equation  is 
modified  to  include  a  dependence  on  loading  direction  and  response  direction  as, 


^  =  wYk^e-,) 


(2.21) 


where  i  is  the  loading  direction  and  j  is  the  response  direction.  The  power  creep  law 
equation  is  modified  as  well  as, 


epjj 


AYf 


(2.22) 


where  i  is  the  loading  direction  and  j  is  the  response  direction. 
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Three  creep  simulations  were  run  at  approximately  75-90%  of  yield  strength  at  a 
constant  load  of  700  MPa  for  10,000  sec  for  the  three  orthogonal  loading  directions. 
The  following  tables  (Tables  2.7,  2.8  and  2.9)  summarize  the  simulation  results  for 
the  three  loading  directions.  Notice  the  anisotropy  between  the  three  cases. 


load  dir. 

i=l 

resp.  dir. 

jj=H 

jj=22 

jj =33 

YS  (MPa) 

917 

n/a 

n/a 

K 

3.9E37 

-9.0E36 

-2.4E35 

n 

18.7 

18.6 

17.8 

A  (1/s) 

2.3E-04 

-8.8E-05 

-1.4E-04 

m 

0.136 

0.145 

0.130 

Table  2.7:  Simulation  and  model  fit  results  showing  anisotropy  (loading  dir.:  i=l). 


load  dir. 

i=2 

resp.  dir. 

jj=H 

jj  =22 

jj =33 

YS  (MPa) 

n/a 

853 

n/a 

K 

-2.2E-27 

4.3E27 

-8.1E25 

n 

14.2 

14.1 

13.3 

A  (1/s) 

-8.4E-05 

1.8E-04 

-1.0E-04 

m 

0.171 

0.205 

0.223 

Table  2.8:  Simulation  and  model  fit  results  showing  anisotropy  (loading  dir.:  i=2). 


The  highly  anisotropic  plastic  response  of  this  material  with  its  dependence  on 
grain  orientation  is  consistent  with  expectations  because  it  is  well  known  that  Ti 
alloys  show  this  property  due  to  the  low  symmetry  of  the  predominant  hep  a  phase 
[42],  The  values  here  for  m  range  between  0.136  to  0.223,  which  is  consistent  with 
experimentally  determined  values  of  approximately  0.2  for  Ti-6A1-4V.  According  to 
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load  dir. 

i=3 

resp.  dir. 

jj=H 

jj  =22 

jj=33 

YS  (MPa) 

n/a 

n/a 

786 

K 

-5.0E16 

-5.4E17 

9.8E17 

n 

9.0 

9.4 

9.4 

A  (1/s) 

-2.5E-04 

-3.0E-04 

5.5E-04 

m 

0.160 

0.178 

0.173 

Table  2.9:  Simulation  and  model  fit  results  showing  anisotropy  (loading  dir.:  i=3). 


Chu  [4],  these  values  for  Ti-6A1-4V  are  significantly  larger  than  those  observed  for 
other  metals  and  alloys  at  low  homologous  temperatures.  This  implies  that  creep 
exhausts  extraordinarily  slowly  in  Ti  alloys  which,  in  general,  is  not  good  for  engi¬ 
neering  applications.  An  understanding  of  how  to  minimize  this  property  in  Ti  alloys 
is  important  and  is  facilitated  by  the  quantitative  understanding  given  in  this  work 
of  how  microstructure  affects  this  ’’cold”  creep  phenomenon. 
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Chapter  3:  Ti-3A1-2.5V  TPS  Material 


The  Space  Shuttles  used  for  NASA’s  recently  retired  Space  Shuttle  program  were 
built  with  a  Thermal  Protection  System  (TPS)  to  withstand  heating  during  atmo¬ 
spheric  reentry.  The  partially  reusable  system  was  composed  mainly  of  ceramic  tiles 
and  blankets  designed  to  be  used  for  multiple  missions,  but  with  a  large  amount  of 
refurbishment  required  between  launches.  In  order  to  drastically  reduce  the  cost  of 
any  future  reusable  launch  vehicle  systems,  metallic  TPS  will  need  to  be  developed. 
There  has  been  a  focus  in  recent  years  for  developing  these  metallic  systems  to  the 
end  of  reducing  the  cost  of  space  access  [3,  2], 

One  challenge  associated  with  developing  these  metallic  TPS  is  accurately  deter¬ 
mining  the  mechanical  properties  of  thin  metallic  foil  materials  in  the  range  of  75-250 
/jrn.  In  [26],  Morrissey  et  al.  describe  a  test  system  that  has  been  developed  that  has 
advanced  the  state-of-the-art  in  testing  capabilities  for  these  thin  sheets  of  material. 
For  the  current  work,  this  test  system  was  used  to  provide  a  large  amount  of  ex¬ 
perimental  data  for  different  specimens  of  thin-gage  rolled  Ti-3A1-2.5V.  Though  the 
most  likely  selection  for  a  titanium  metallic  TPS  alloy  is  a  high  temperature  capable 
material  such  as  Ti-6Al-2Mo-4Zr-2Sn,  the  readily  available  Ti-3A1-2.5V  commercial 
alloy  was  selected  as  a  suitable  substitute  material  to  demonstrate  the  effect  of  cold 
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rolling  on  the  texture  of  the  processed  material  and  the  corresponding  anisotropic 
mechanical  response  of  the  sheet. 

Following  a  similar  process  described  in  Chapter  2,  the  following  chapter  describes 
the  modeling  work  completed  using  the  mechanical  response  experimental  data  as 
well  as  the  micrographs  and  EBSD  scans  of  the  pre-tested  material.  Because  of  the 
pronounced  effect  of  the  rolling  on  the  material’s  microstructure  and  the  large  amount 
of  test  data  for  specimens  of  different  orientations,  the  calibration  and  validation  of 
the  highly  anisotropic  material  parameters  for  this  material  were  able  to  be  completed 
with  a  good  degree  of  fidelity. 

3.1  Experimental  Set-up  and  Results 

In  comparison  with  the  experimental  data  available  for  the  Ti-6A1-4V  material 
used  in  Chapter  2,  for  this  material  an  abundance  of  experimental  data  was  available 
for  model  comparison,  calibration,  and  validation.  Using  a  newly  developed  test 
system  for  mechanical  characterization  of  thin  metallic  foils  discussed  by  Morrissey 
et  al.  in  [26],  mechanical  response  data  for  a  number  of  foil  specimens  was  provided. 
See  Figure  3.5  for  a  picture  of  the  experimental  set-up. 

Foil  sample  specimens  were  taken  from  the  bulk  sheet  material  cut  to  be  aligned 
with  both  the  rolling  direction  (RD  specimens)  and  the  transverse  direction  (TD 
specimens).  See  Figure  3.1  for  a  schematic.  On  these  individual  specimens,  strain 
gages  were  attached  to  collect  local  strain  data  in  both  the  axial  direction  (0°  strain 
gages)  and  transverse  direction  (90°  strain  gages).  Some  tests  included  45°  strain 
gages.  See  Figure  3.2  for  a  schematic  showing  the  strain  gage  placements.  A  schematic 
of  a  typical  foil  specimen  with  dimensions  is  shown  to  scale  in  Figure  3.3.  Thicknesses 
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of  the  foil  specimens  are  4.5  MIL  (approximately  120  /mi).  See  Figure  3.4  for  a  close- 
up  picture  of  one  of  the  experimental  specimens  showing  strain  gage  placement  for 
the  0°  and  90°  strain  gages.  Notice  wires  running  to  the  back  side  of  the  specimen;  the 
reverse  side  of  the  specimens  have  a  similar  strain  gage  set-up.  This  was  on  account 
of  the  thin  nature  of  the  specimens  in  order  to  confirm  that  applied  displacements 
and  loads  were  not  causing  any  warping  of  the  specimen  in  out-of-plane  directions. 
Similar  front  and  back  strain  gage  measurements  confirmed  a  well-executed  test. 


Figure  3.1:  Schematic  of  the  bulk  sheet  material  showing  how  specimens  are  cut  at 
different  orientations.  Not  to  scale. 

Various  specimens  were  tested  under  different  loading  conditions.  Tensile  results 
were  obtained  by  displacement  controlled  tests  for  both  RD  specimens  and  TD  spec¬ 
imens.  In  addition,  specimens  from  both  directions  were  tested  for  creep  and  dwell 
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< -  loading  direction  - ► 


Figure  3.2:  Schematic  of  a  specimen  showing  strain  gage  placement  and  defining  the 
specimen  coordinate  system  (1-2)  with  1  being  the  rolling  direction.  Not  to  scale. 
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Figure  3.3:  Specimen  schematic  to  scale  with  dimensions  in  mm. 
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Figure  3.4:  Close-up  picture  of  a  specimen  fitted  in  rig  for  testing  at  AFRL/RX. 
Strain  gages  for  this  specimen  are  in  the  0°  and  90°  orientations. 


Figure  3.5:  Picture  of  experimental  set-up  at  AFRL/RX. 
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fatigue  results  through  load  controlled  tests.  The  load  controlled  tests  were  conducted 
after  the  displacement  controlled  tests  helped  to  determine  yield  strength.  A  number 
of  loads  were  tested  around  the  range  of  95%  of  yield.  Constant  creep  load  and  dwell 
loading  pattern  tests  were  conducted  with  the  same  maximum  load  for  direct  com¬ 
parison  of  the  dwell  fatigue  effect  on  mechanical  response.  See  Table  3.1  for  a  detailed 
summary  of  the  various  test  specimens.  See  Figure  3.6  for  experimental  tensile  results 
and  Figure  3.8  for  experimental  creep  and  dwell  results.  For  the  dwell  experiments, 
strain  values  are  reported  at  the  midpoint  during  the  dwell  duration.  See  Figure  3.7 
for  a  plot  of  the  load  pattern  for  1  dwell  cycle. 

3.2  CPFEM  Simulation  of  Polycrystalline  Ti-3A1-2.5V 

In  order  to  accurately  simulate  the  micromechanical  response  of  polycrystalline 
materials,  it  is  important  to  develop  both  a  good  microstructure  model  as  well  as  a 
good  micromechanical  model.  These  separate  but  related  aspects  to  the  simulation 
are  discussed  in  the  following  pages  beginning  with  the  microstructural  simulation 
model. 

3.2.1  3D  Polycrystalline  Microstructural  Simulation  and  Mesh 
Generation 

The  following  section  discusses  the  2D  characterization  and  generation  of  3D  mi¬ 
crostructure  models  from  the  2D  data.  Section  2.1.1  gave  a  fuller  set  of  details  for 
the  procedure  while  this  particular  section  will  focus  on  the  extension  of  that  pro¬ 
cedure  to  include  using  EBSD  scans  from  3  orthogonal  surfaces  for  a  more  precise 
representation  of  the  morphological  anisotropy  of  this  rolled  material. 
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Figure  3.6:  Experimental  tensile  res 
0°  strain  gage  corresponds  to  i=j=l, 
strain  gage  corresponds  to  i=l,  j=2. 


Its  (full  results  (a)  and  a  close-up  (b))  where 
)0°  strain  gage  corresponds  to  i=j =2,  and  45° 
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Strain  values  in  the  strain 
versus  time  plots  correspond 


Figure  3.7:  Dwell  loading  pattern  for  1  cycle  showing  midpoint  where  data  was  col¬ 
lected  past  10  cycles.  Figure  generated  by  AFRL/RX  with  some  modification.  *See 
Table  3.1  for  specific  ramp  rates. 


Microstructural  Characterization  of  the  a  Ti-3A1-2.5V 

Three  EBSD  scans  of  the  Ti-3A1-2.5V  material  taken  by  AFRL/RX  from  3  or¬ 
thogonal  surfaces  were  used  to  characterize  the  microstructure.  A  schematic  of  the 
rolled  material  labeling  the  three  orthogonal  surfaces  was  given  in  Figure  3.1  and  the 
three  EBSD  scans  are  shown  in  Figure  3.9. 

In  contrast  to  the  bi-modal  a-f3  Ti-6A1-4V  material  studied  in  Chapter  2,  this 
polycrystalline  material  is  only  composed  of  a  phase  which  has  an  hep  crystal  struc¬ 
ture  and  associated  slip  systems  (see  Figure  2.2).  Though  this  removes  some  of  the 
modeling  complexity  (not  needing  to  model  the  complicated  transformed  ft  phase), 
the  rolled  nature  of  this  material  gives  rise  to  a  pronounced  crystallographic  texturing 
as  well  as  elongated  grain  shapes  which,  as  will  be  shown,  have  a  significant  influence 
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Figure  3.8:  Experimental  creep  and  dwell  results  for  (a)  0°  strain  gages  and  (b)  90° 
strain  gages. 
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Figure  3.9:  EBSD  scans  of  (a)  surface  face,  (b)  transverse  face,  and  (c)  longitudinal 
faces;  with  (d)  color  legend. 
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on  mechanical  response  properties.  It  is  important  that  these  features  are  accurately 
represented  by  the  computational  models. 

Using  the  code  described  in  [8],  the  EBSD  scans  from  the  three  orthogonal  faces 
were  analyzed  to  determine  a  number  of  crystallographic  and  morphological  mi- 
crostructural  characteristics  including  crystal  size  distribution,  aspect  ratio  distribu¬ 
tion  and  principal  axis  orientation  distribution  (assuming  ellipsoidal-shaped  grains), 
grain  size  versus  number  of  neighbor  distribution,  crystallographic  orientation  distri¬ 
bution,  misorientation  distribution,  and  microtexture  distribution. 

Microstructural  Simulation  Procedure 

Using  a  similar  microstructural  simulation  procedure  as  described  in  Chapter  2, 
a  representative  synthetic  microstructure  is  generated  that  is  statistically  equivalent 
to  the  real  microstructure  (as  characterized  by  the  statistics  given  in  the  previous 
section).  As  mentioned,  the  rolled  nature  of  this  material  gives  rise  to  deformed  grain 
shapes.  The  rolled  surface  scan  and  aspect  ratio  analysis  shows  grains  that  have  been 
”  flattened”  while  the  longitudinal  face  shows  grains  that  have  been  elongated  in  a 
”  bread-stick”  fashion.  One  significant  difference  between  the  modeling  work  for  this 
rolled  material  as  compared  to  that  for  the  Ti-6A1-4V  material  is  that  morphological 
aspect  ratio  data  from  the  three  orthogonal  faces  was  matched  in  the  full  3D  synthetic 
microstructure.  Grain  shapes  assumed  to  be  elliptic  in  2D  are  further  assumed  to  be 
ellipsoidal  in  3D  with  principal  axis  length  distributions  individually  matched  for 
each  2D  face.  Overlapping  of  grains  and  empty  space  is  assigned  to  the  nearest 
grains  by  an  iterative  process  where  statistics  are  monitored  until  convergence  is 
reached  and  the  entire  cubic  volume  is  filled  with  no  overlapping.  After  the  grain 
shapes  have  been  placed  and  the  volume  has  been  filled,  orientations  are  assigned  to 
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each  grain  based  off  a  random  sampling  of  the  orientation  distribution  function  and 
orientations  are  swapped  until  convergence  is  reached  through  an  iterative  process 
where  misorientation  and  microtexture  distribution  errors  are  monitored. 

A  synthetic  microstructure  comprised  of  534  grains  was  built  in  the  fashion  de¬ 
scribed  above  and  can  be  considered  statistically  equivalent  with  respect  to  these 
important  statistics.  Chapter  2  has  shown  that  the  microstructure  simulation  al¬ 
gorithm  converges  with  respect  to  grain  size  distribution  for  structures  larger  than 
approximately  300  grains.  Figures  3.10  and  3.11  show  the  534-grain  structure  from 
4  viewpoints.  Notice  the  difference  in  grain  shapes  from  the  3  different  faces  of  the 
microstructure.  Visual  inspection  shows  that  the  synthetic  generator  algorithm  has 
faithfully  reproduced  a  rolled  microstructure. 


ND 


RD 


TD 


Figure  3.10:  Voxclized  volume  for  the  534-grain  microstructure.  The  cube  length 
dimension,  /0,  is  24.6  /irn. 


Grain  size  distribution  is  another  important  microstructural  feature  that  has  a 
strong  effect  on  mechanical  properties  especially  as  defined  by  the  Hall-Petch  rela¬ 
tionship  2.14.  Figure  3.12  shows  the  grain  size  distribution  comparison  for  the  real 
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(a) 


(b) 


(c) 


Figure  3.11:  2D  views  of  the  synthetic  structure’s  (a)  surface  face,  (b)  transverse  face, 
and  (c)  longitudinal  face.  Notice  ’’flattened”  shapes  on  surface  face  and  elongated 
shapes  on  longitudinal  face. 


microstructure  as  compared  to  the  synthetic  microstructure.  The  version  of  the  syn¬ 
thetic  generator  code  available  at  the  time  had  hard-coded  standard  deviation  cut-offs 
that  limited  the  maximum  and  minimum  grain  sizes  that  could  be  sampled  from  the 
lognormal  curve  fit.  The  grain  size  distribution  based  on  volume  fraction  (Figure  3.12 
(b))  helps  visualize  the  pronounced  effect  of  not  matching  the  upper  tail.  The  reason 
for  the  discrepancy  is  that  the  synthetic  structure  does  not  have  as  many  of  the  largest 
sized  grains  as  the  real  microstructure  has.  See  Section  3.3.1  for  a  short  discussion  on 
a  possible  implication  of  this  discrepancy.  It  is  recommended  that  future  work  should 
investigate  through  a  sensitivity  the  effect  of  not  matching  the  upper  tail  of  the  grain 
size  distribution. 

Another  very  important  microstructural  statistic  that  effects  mechanical  response, 
is  crystallographic  orientation  distribution.  This  characteristic  is  particularly  impor¬ 
tant  for  the  current  material  because  of  the  anisotropic  nature  of  hep  crystals  (which 
comprise  100%  of  the  structure)  and  the  highly  textured  nature  of  the  rolled  material. 
Figure  3.13  gives  a  comparison  of  a  (0002)  pole  figure  for  the  sample  data  acquired 
from  the  EBSD  scans  and  the  534  grains  of  the  synthetic  structure.  A  pole  figure 
point  density  (PFPD)  distribution  comparison  is  included  for  a  more  quantitative 
comparison  of  the  real  structure’s  crystal  orientation  and  the  synthetic  structure’s 
crystal  orientations.  The  crystallographic  texture  of  the  rolled  Ti-3A1-2.5V  is  very 
similar  to  that  seen  for  rolled  commercially  pure  titanium  as  reported  by  Lutjering 
and  Williams  [25].  According  to  Codings  [5],  it  is  agreed  that  such  a  texture  results 
from  a  competition  between  slip  which  rotates  the  basal  poles  toward  the  sheet  nor¬ 
mal,  and  twinning  which  tends  to  rotate  them  into  the  transverse  direction.  Hence, 
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(a) 


(b) 


Figure  3.12:  Grain  size,  D,  distribution  comparison  based  on  (a)  number  fraction  of 
grains  (IV/)  and  (b)  volume  fraction  of  grains  (V/). 
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the  (0002)  pole  figure  for  rolled  titanium  shows  a  concentration  along  the  TD  axis 
with  some  preference  to  the  ND  origin. 


(a)  (b) 


Figure  3.13:  a  phase  (0002)  pole  figures  for  (a)  Sample  data  and  (b)  Synthetic  struc¬ 
ture  as  well  as  (c)  pole  figure  point  density  (PFPD)  distributions  of  the  2  pole  figures 
for  a  more  quantitative  comparison  of  the  crystallographic  matching. 


The  synthetic  generator  code  also  matches  the  more  local  texturing  metrics  of 
misorientation  distribution  and  microtexture  distribution.  See  Figures  3.14  and  3.15 
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for  a  comparison  of  these  statistics  for  the  real  structure  and  synthetic  structure.  Al¬ 
though  not  as  important  for  this  particular  study,  matching  of  these  metrics  becomes 
especially  important  for  more  local  mechanical  responses  such  as  crack  nucleation 
which  is  part  of  the  proposed  future  work  (see  Chapter  4). 


Figure  3.14:  Misorientation  distribution  (MoDF)  comparison.  GBf  corresponds  to 
either  grain  boundary  length  fraction  or  grain  boundary  surface  area  fraction. 


Mesh  Generation 

The  commercial  software  package  Simmetrix  [37]  was  used  to  mesh  the  voxelized 
volume  generated  by  the  synthetic  structure  generator  code.  The  meshed  structure 
consists  of  151,474  four-noded  tetrahedral  elements  and  a  total  of  26,437  nodes.  See 
Figure  3.16  (b)  for  a  3D  view  of  the  meshed  structured.  The  mesh  is  checked  for 
distorted  elements  and  it  is  seen  that  less  than  0.01%  of  the  elements  have  an  aspect 
ratio  of  40  or  higher.  Section  2.1.1  has  shown  that  for  the  Ti-6A1-4V  material  meshes 
of  500  or  600  grains  that  are  meshed  with  similar  resolution  do  not  show  any  significant 
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Figure  3.15:  Microtexture  distribution  comparison  where  microtexture  is  defined  by 
the  number  fraction  of  neighbors  ( Nf n)  with  misorientation  less  than  15°.  Nf  is  the 
number  fraction  of  grains. 

mesh  dependence  when  it  comes  to  the  volume  averaged  creep  response.  Because  the 
current  work  is  primarily  concerned  with  volume  averaged  macroscopic  responses, 
it  is  determined  that  the  overall  response  will  not  be  greatly  affected  by  any  mesh 
dependencies  for  this  534-grain  structure. 

This  microstructure  and  mesh  will  be  referred  to  in  this  work  as  Microstructure 
A  or  MS -A  for  short. 

Simplified  Microstructure  for  Model  Updates,  Sensitivity  Studies  and  Pa¬ 
rameter  Calibration 

Because  of  the  large  amount  of  computational  expense  of  the  meshed  534-grain 
synthetic  structure  (MS- A),  a  second  simplified  microstructure  and  mesh  were  de¬ 
veloped  with  some  justifiable  simplifications  for  the  purpose  of  performing  model 
updates,  sensitivity  studies  and  parameter  calibration. 


63 


Figure  3.16:  Voxelized  volume  (a)  before  and  (b)  after  meshing  for  MS- A.  The  cube 
length  dimension,  1$,  is  24.6  /ini. 

A  regular  grid  structure  consisting  of  2744  cubic  grains  with  8  tetrahedral  el¬ 
ements  per  grain  (for  a  total  of  16464  elements)  and  3375  nodes  was  generated. 
Orientations  were  assigned  to  the  grains  from  a  random  sampling  of  orientations 
from  the  EBSD  scan  of  the  real  microstructure  producing  a  microstructure  that  was 
statistically  equivalent  in  terms  of  the  orientation  distribution  function. 

According  to  an  unpublished  work  by  Anahid  and  Ghosh  on  a  macroscopic  yield 
surface  model  for  titanium,  it  is  not  as  important  for  the  current  work  to  consider 
the  misorientation  and  microtexturing  distributions.  The  effect  of  these  character¬ 
istics  on  the  macroscopic  volume  averaged  constant  strain  rate  and  creep  responses 
are  not  as  pronounced  as  the  relatively  larger  effects  of  grain  size  and  orientation 
distribution.  Again,  as  mentioned  above,  local  effects  that  are  highly  dependent  on 
local  neighborhoods  of  grains  are  more  sensitive  to  metrics  such  as  misorientation 
and  microtexturing.  Crack  nucleation  models,  for  example,  require  these  characteris¬ 
tics  of  the  real  sample  to  be  accurately  reproduced  in  any  synthetic  representations. 
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Since  the  scope  of  the  current  work  is  not  focused  on  investigating  local  phenomena, 
matching  of  misorientation  and  texture  distribution  was  neglected  in  this  simplified 
microstructure  for  the  model  updates,  parameter  sensitivity  study  and  calibration 
efforts. 

Even  though  the  finite  element  mesh  has  uniform  grain  sizes  (1  /ini  cubes),  the 
Hall-Petch  size-effect  as  well  as  a  newly  developed  aspect  ratio  dependence  (see  Sec¬ 
tion  3.2.2)  is  incorporated  into  the  model  by  assigning  different  D  values  to  the  grains 
based  on  grain  size  and  aspect  ratio  distributions  of  the  real  specimens.  So,  although 
the  FEM  mesh  has  grains  of  uniform  size  and  shape,  individual  grains  plastically  de¬ 
form  differently  from  each  other  based  on  their  D  value  assignment  which  corresponds 
directly  to  the  real  material. 

In  the  simplified  cubic  grain  structure  for  the  H-aspect  ratio  dependence  (discussed 
in  Section  3.2.2),  grains  are  oriented  with  their  principal  axes  aligned  with  the  rolling 
direction.  This  assumption  is  reasonable  considering  the  ellipsoidal  principal  axes  of 
the  grains  of  the  rolled  structure  are  almost  all  aligned  with  the  rolling  direction  with 
only  slight  variation  as  seen  by  visual  inspection  of  the  longitudinal  face  (Figure  3.11 

(c))- 

This  microstructure  and  mesh  will  be  referred  to  in  this  paper  as  Microstructure 
B  or  MS-B  for  short. 

Table  3.2  summarizes  the  differences  between  the  real  microstructure  and  the 
synthetic  structures  MS-A  and  MS-B. 
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Figure  3.17:  Regular  grid  model  of  MS-B  ((a)  before  and  (b)  after  meshing).  The 
cube  length  dimension,  /0,  is  14.0  prn. 


real 

MS-A 

MS-B 

av.  grain  size  (no.  fract.),  D 

3.07  pm 

3.08  pm 

3.17  pm 

st.  dev.  grain  size  (no.  fract.) 

1.77  pm 

1.52  pm 

1.50  pm 

av.  grain  size  (vol.  fract.) 

7.08  pm 

5.19  pm 

5.22  pm 

st.  dev.  grain  size  (vol.  fract.) 

3.20  pm 

1.52  pm 

1.54  pm 

max.  grain  size 

16.8  pm 

7.35  pm 

7.37  pm 

min.  grain  size 

1.42  pm 

0.71  pm 

0.74  pm 

vol.  fract.  ol  phase 

100% 

100% 

100% 

ellipsoid  princ.  axes  matched? 

n/a 

yes 

roll.  dir. 

orientation  dist.  matched? 

n/a 

yes 

yes 

misorientation  dist.  matched? 

n/a 

yes 

no 

texture  dist.  matched? 

n/a 

yes 

no 

mesh  approx. 

n/a 

real,  shapes 

cubes 

plastic  property  D  approx. 

n/a 

equi.  spheres  equi.  ellip. 

no.  of  grains 

—2  billion 

534 

2,744 

no.  of  elements 

n/a 

151,474 

16,464 

no.  of  elements  per  grain 

n/a 

-284 

8 

no.  of  nodes 

n/a 

26,437 

3,375 

Table  3.2:  Summary  of  real  microstructure  compared  with  MS-A  and  MS-B. 
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3.2.2  Micromechanical  Analysis  of  Simulated  Polycrystalline 
Microstructures  Using  a  Crystal  Plasticity  Constitu¬ 
tive  Model 

The  isothermal,  size-dependent  and  rate-dependent  crystal  plasticity  finite-element 
model  discussed  in  Section  2.1.2  is  used  along  with  an  in-house  parallel  code  to  sim¬ 
ulate  the  mechanical  responses  of  MS-A  and  MS-B.  The  initial  elastic  and  plastic 
parameters  used  were  calibrated  by  previous  single  crystal  work  discussed  in  [13]  and 
[7]  where  a  Ti-6A1  single  crystal  experimental  data  was  available  for  calibration  of 
individual  slip  system  parameters.  See  Section  2.1.2  for  initial  values  and  Section 
3.3.1  for  the  calibration  efforts  for  this  particular  alloy.  An  implicit  backward  Euler 
time-integration  scheme  is  used  for  the  solution  to  the  dynamic  problem.  An  update 
to  the  CPFEM  model  is  discussed  and  proposed  in  the  following  section. 

Crystal  Plasticity  Constitutive  Model  Update  to  Include  Grain  Aspect 
Ratio  Dependent  Yielding 

The  constitutive  equations  governing  the  material  behavior  are  all  given  in  Section 
2.1.2.  The  one  notable  difference  being  that  the  Ti-3A1-2.5V  material  has  no  trans¬ 
formed  j3  colonies  and  therefore  no  bcc  phase  to  model.  This  reduces  the  complexity 
of  the  constitutive  equations;  however,  the  rolled  nature  of  this  material  gives  rise  to 
the  need  for  an  update  to  the  model  to  account  for  grain  aspect  ratio  dependency  on 
the  size  effect  for  the  hardness  law,  given  previously  as  Equation  2.11.  As  mentioned, 
the  slip  system  deformation  resistance,  g°  for  each  slip  system  is  given  a  Hall-Petch- 
type  size  effect  as  ga  =  +  ^=.  In  previous  work  (including  the  work  done  for 

the  Ti-6A1-4V  material  in  Chapter  2),  the  value  of  the  characteristic  length  scale,  Da 
assumed  that  each  grain  was  spherical  in  shape.  Therefore,  there  was  no  aspect  ratio 
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dependence  governing  the  size  effect  on  the  ga.  This  is  a  reasonable  assumption  as 
long  as  grain  and  colony  shapes  are  relatively  spherical  in  shape.  However,  for  rolled 
materials,  this  assumption  is  no  longer  reasonable. 

An  update  for  the  model  is  proposed  and  implemented  which  assigns  different 
Da  values  for  each  slip  system  (30  per  grain  for  hep  crystals)  based  on  the  distance 
from  the  origin  of  an  equivalent  ellipsoid  to  a  point  on  the  ellipsoid’s  surface  in  the 
direction  of  slip  multiplied  by  2.  A  schematic  of  an  ellipsoid  is  shown  in  Figure 
3.18  with  the  new  characteristic  length  given  as  D.  This  update  to  the  model  gives 
a  more  physical  representation  of  the  size  effect  on  the  yield  property  of  individual 
grains  which  was  originally  developed  based  on  the  fact  that  grain  boundaries  impede 
dislocation  motion. 


Figure  3.18:  Schematic  of  an  equivalent  ellipsoid  representing  a  grain  showing  how 
the  new  characteristic  length  scale,  Da,  is  constructed  based  on  the  intersection  of 
a  line  from  the  origin  in  the  slip  direction  to  the  intersection  of  the  surface  of  the 
ellipsoid. 


In  actuality,  for  moderately  rolled  sheets  of  polycrystalline  material  where  the 
grains  are  seen  to  be  elongated  by  approximately  1.5  times  (as  is  the  case  with  the 
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current  material),  it  is  seen  that  the  effect  of  adding  this  dependence  is  not  signif¬ 
icant  in  the  macroscopic  volume  averaged  constant  strain  rate  and  creep  response 
parameters.  It  is  conjectured  that  this  functionality  would  be  very  important  for  the 
more  locally  sensitive  crack  nucleation  response  of  the  material  (see  Chapter  4  for 
a  discussion  on  future  work).  In  order  to  see  the  effect  of  the  model  update  in  the 
macroscopic  response,  constant  strain  rate  simulations  were  run  on  a  microstructure 
where  grains  were  elongated  with  an  average  a/b  aspect  ratio  of  4  and  a/c  aspect  ratio 
of  16,  where  a,  b,  and  c  correspond  to  the  principal  axis  lengths.  The  microstructure 
was  loaded  in  the  rolling  direction  (the  direction  of  the  elongated  grains).  Comparison 
of  the  simulation  results  before  and  after  the  model  update  gave  a  3.6%  difference  in 
yield  strength  (as  measured  by  the  0.2%  strain  offset  method).  The  aspect  ratio  de¬ 
pendent  results  predict  a  higher  yield  strength  than  that  of  the  model  which  assumes 
equivalent  spheres. 

Boundary  and  Loading  Conditions 

Because  of  the  limitation  of  micromechanical  CPFEM  to  analyze  microstructures 
with  a  maximum  number  of  grains  of  the  order  of  approximately  300-3000  grains,  the 
absolute  dimensions  of  the  simulation  models  need  to  be  limited  to  the  order  of  about 
50-150yum.  This  is  much  smaller  than  the  actual  experimental  specimens  which  have 
gage  lengths  containing  approximately  2  billion  grains  (as  shown  in  Table  3.2)  and 
dimensions  reaching  the  order  of  centimeters  (see  Figure  3.3).  Boundary  conditions 
are  used  to  approximate  the  behavior  of  the  boundaries  of  these  micro- cubes. 

For  the  constant  strain  rate  simulations,  in  order  to  suppress  rigid  body  modes, 
symmetric  constraint  conditions  are  applied  where  the  back  faces  of  the  cube  are 
constrained  as  u*=0  on  the  ith-face,  where  i  goes  from  1  to  3.  Directions  1  and  2 
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represent  specimen  orientations  in  the  RD  and  TD  orientations  while  direction  3  is 
the  sheet  normal  direction.  The  cube  is  free  to  move  in  all  other  directions.  One 
of  the  outer  faces  is  given  a  constant  strain  rate  displacement  boundary  condition 
of  Unit)  =  ^(exp^jt)  —  1),  where  l0  is  the  initial  dimension  of  the  cube  and  e-  is 
the  applied  constant  strain  rate  which  for  these  experiments  is  taken  to  be  9.76E- 
5  1/s.  For  the  actual  experiments,  a  constant  stroke  displacement  rate  of  6.16E-3 
nnn/s  gave  rise  to  an  average  strain  rate  of  7.63E-5  1/s  for  the  elastic  portion  of 
the  experiment  and  an  average  strain  rate  of  9.76E-5  1/s  for  the  plastic  portion  of 
the  experiment  (the  portion  following  yield).  Different  simulations  were  run  using 
the  average  elastic  rate  for  the  entire  simulation  and  the  average  plastic  rate  for 
the  entire  simulation  and  very  little  difference  was  seen  in  the  macroscopic  response; 
therefore,  for  simplification,  the  assumption  of  9.76E-5  1/s  as  the  strain  rate  for  the 
entire  simulation  is  taken  to  be  reasonable.  It  is  conjectured  that  for  higher  strain 
rate  experiments  that  this  assumption  would  no  longer  be  valid  and  a  more  precise 
representation  of  the  applied  strain  rates  would  be  necessary  to  accurately  model 
the  experimental  results.  Constant  strain  rate  simulations  are  run  to  about  5%  total 
strain  to  ensure  that  the  full  tensile  experiment  range  is  simulated. 

For  the  creep  simulations,  both  symmetric  boundary  conditions  (as  discussed 
above)  and  minimum  boundary  conditions  were  used  to  ensure  that  boundary  con¬ 
ditions  were  not  affecting  the  macroscopic  response  too  much.  Minimum  boundary 
conditions  correspond  to  the  back  face  of  the  cube  being  constrained  from  motion 
in  the  loading  direction,  1  line  on  this  same  back  face  in  the  normal  direction  being 
constrained  from  motion  in  the  non-loading  direction  (to  constrain  rotational  motion 
about  the  loading  direction),  and  1  point  at  the  origin  being  constrained  in  all  three 
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directions  to  restrain  the  rigid  motion  in  the  normal  direction.  It  was  seen  for  MS-B 
that  using  symmetric  or  minimum  boundary  conditions  for  the  creep  simulations  did 
not  affect  the  macroscopic  response  significantly.  This  verifies  that  we  are  not  hav¬ 
ing  a  boundary  condition  dependence  for  these  runs.  For  creep,  constant  loads  are 
applied  to  one  of  the  outer  faces  for  about  334  hours. 

3.3  CPFEM  Simulation  Results  with  Material  Parameter 
Calibration  Procedure  and  Model  Validation  Using  Ex¬ 
perimental  Results 

Using  the  material  parameters  calibrated  for  previous  titanium  alloys,  initial  ten¬ 
sile  simulation  results  for  the  Ti-3A1-2.5V  material  did  not  match  well.  This  is  under¬ 
standable  because  this  is  a  different  alloy  with  a  different  processing  history.  Yield 
strengths  near  900  MPa  were  predicted  for  specimens  of  both  RD  and  TD  orientations 
while  experimental  results  are  nearer  600  MPa.  In  addition,  the  post-yield  slopes  were 
not  matching  well  and  the  elastic  portion  of  the  simulations  were  also  off  by  about 
10%.  Fortunately,  the  large  amount  of  experimental  data  available  allowed  for  some 
of  the  data  to  be  used  for  calibration  purposes. 

The  displacement  controlled  tensile  experiments  are  used  for  calibration.  Once 
the  calibration  is  completed  for  the  constant  strain  rate  simulations  using  the  tensile 
experimental  results,  validation  of  the  model  is  achieved  by  predicting  creep  results 
using  the  calibrated  model. 

3.3.1  Material  Parameter  Calibration  Using  Displacement 
Controlled  Experimental  Data  and  Sensitivity  Studies 

Because  of  the  anisotropic  nature  of  the  clastic  and  crystal  plasticity  constants, 
calibration  efforts  required  the  use  of  experimental  data  for  both  the  RD  and  TD 
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specimens  to  fully  calibrate  these  directionally  dependent  constants.  The  hep  crystals 
exhibit  different  responses  in  the  different  directions  and  it  was  seen  that  the  RD  and 
TD  specimens  were  sensitive  to  different  material  constants  to  different  degrees.  For 
visualization  of  slip  family  activity,  Figure  3.19  gives  the  distributions  of  softness 
indices,  SI,  as  defined  by  Equation  2.20  where  here  softness  indices  are  shown  for  the 
specific  slip  families  rather  than  for  the  softest  slip  system  of  the  grain.  ’’Activity” 
corresponds  to  values  of  softness  index  above  1  (the  black  vertical  line  in  Figure 
3.19).  The  simulations  shows  larger  basal  slip  system  activity  in  the  TD  specimens 
(red,  dotted  curve)  than  in  the  RD  specimens  where  there  is  almost  no  basal  activity 
(red,  solid  curve).  In  addition,  the  RD  specimens  (green,  solid  curve)  show  more 
prismatic  activity  than  do  the  TD  specimens  (green,  dotted  curve).  The  difference 
in  the  plastic  response  of  these  slip  systems  is  what  largely  accounts  for  the  different 
macroscopic  responses  in  the  polycrystalline  material.  For  this  reason,  it  would  not 
be  possible  to  calibrate  using  RD  specimen  data  and  then  predict  the  TD  specimen 
response  accurately. 

Elastic  Parameter  Calibration 

The  elastic  parameters  were  calibrated  through  sensitivity  studies  where  the  5 
independent  constants  were  varied  while  holding  all  others  constant.  Using  the  sen¬ 
sitivity  study  results,  the  parameters  were  tweaked  until  the  elastic  portion  of  the 
constant  strain  rates  simulations  for  both  RD  and  TD  specimens  matched  tensile 
experimental  results.  Table  3.3  gives  the  final  calibrated  results  and  Table  3.4  shows 
a  comparison  of  the  experimental  elastic  properties  versus  the  simulated  properties 
(before  and  after  calibration). 
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—  -TDspec.;  basal 
- RD  spec;  basal 

—  TD  spec.:  prism 
- RD  spec.;  prism 

—  -TDspec.;  pyr.  <a> 
- RD  spec.;  pyr.  <a> 

—  -TDspec.;  pyr.  <c+a> 
- RD  spec.;  pyr.  <c+a> 


Figure  3.19:  Distribution  of  softness  index  (SI  =  SFa  *  YS/g q  )  of  slip  families.  Nf 
is  the  number  fraction  of  grains.  Values  of  SI  above  1  correspond  to  slip  activity. 


Cij  parameter 

Value  (GPa) 

C11—C22 

162.2 

C33 

95.6 

C\2=C'2l 

83.9 

C13  =Czi  =<7-23  =  C*32 

204.0 

C44 

Cu-Cn 

C55=Ce6 

87.0 

Other 

0 

Table  3.3:  Calibration  of  the  Ti-3A1-2.5V  transversely  isotropic  elastic  parameter 
tensor. 


Exp. 

Sim. 

Sim. 

spec. 

Pre-calib. 

Post-calib. 

orient. 

Eit  (GPa) 

Vn 

Eu  (GPci)  Va 

Ea  (GPr)  Vn 

RD  (i=l) 

102 

0.390 

109 

0.376 

102 

0.383 

TD  (i=2) 

111 

0.417 

121 

0.416 

110 

0.412 

Table  3.4:  Elastic  response  parameter  comparison  for  experiments  with  the  simula¬ 
tions  before  and  after  calibration. 
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Crystal  Plasticity  Parameter  Calibration 


Crystal  plasticity  parameter  calibration  was  accomplished  by  the  aid  of  a  wide 
range  of  sensitivity  studies  performed  on  the  MS-B  microstructure.  Studied  parame¬ 
ters  include:  h0,  7,  m,  qab ,  n,  r,  and  initial  slip  system  hardness  parameters,  ,  for  slip 
in  the  basal  <a>,  prismatic  <a>,  pyramidal  <c>,  and  pyramidal  <c+a>directions. 
Additionally,  a  small  model  update  was  made  where  the  g  saturation  hardness  pa¬ 
rameter  is  defined  as  a  scaled  up  value  from  the  size-effect-modified  go  values  given  as 
go-  This  update  is  a  safe  guard  that  the  size-effect  modifications  on  the  g$  values  do 
not  create  an  unphysical  situation  where  the  saturation  value  is  actually  less  than  the 
initial  hardness  of  the  slip  system.  This  was  seen  to  occur  with  the  previous  model 
causing  an  unphysical  softening  of  slip  systems  as  plastic  deformation  increased. 

Table  3.5  gives  the  previously  calibrated  values  as  well  as  the  calibrated  values 
for  MS-A.  When  the  calibrated  values  for  MS-A  were  applied  to  MS-B ,  the  results 
did  not  match  as  well  as  would  be  desired.  The  reason  for  this  can  be  explained 
by  the  differences  between  these  2  microstructures  as  outlined  in  Table  3.2.  A  small 
amount  of  additional  calibration  was  required  to  match  the  simulation  results  for  MS- 
B  (as  shown  in  the  last  column  of  3.5).  This  issue  was  not  seen  in  an  unpublished 
work  by  Anahid  and  Ghosh  based  on  a  comparison  of  a  cube-shaped  grain  model 
and  a  real-shaped  grain  model  of  the  Ti-6A1-4V  alloy  discussed  in  Chapter  2.  It  is 
conjectured  that  modeling  of  the  actual  grain  shapes  has  a  pronounced  effect  on  the 
CPFEM  results  in  this  case  because  of  the  elongated  grain  shapes  for  the  current 
rolled  material. 

Figure  3.20  shows  the  MS-A  constant  strain  rate  simulation  results  after  calibra¬ 
tion  as  compared  with  experimental  results  and  Figure  3.21  gives  the  results  for  MS-B 
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Prev. 

Ti-3A1-2.5V 

Parameter 

Unit 

Calib. 

MS-A 

MS-B 

go  basal  <a> 

(MPa) 

385.0 

239.0 

199.2 

go  prism  <a> 

(MPa) 

365.0 

175.6 

146.4 

g0  pyr.  <a> 

(MPa) 

n/a 

289.1 

240.9 

go  pyr.  <c+a> 

(MPa) 

623.3 

456.1 

380.1 

g  basal  <a> 

(MPa) 

450.0 

<70*1-2 

<70*1-2 

g  prism  <a> 

(MPa) 

550.0 

<70*1-6 

<70*1-6 

9  Pyr-  <a> 

(MPa) 

1650.0 

<70*1-6 

<70*1-6 

g  pyr.  <c+a> 

(MPa) 

450.0 

<70*1-6 

<70*1-6 

ho  basal  <a> 

100.0 

100.0 

100.0 

ho  prism  <a> 

100.0 

100.0 

100.0 

ho  pyr.  <a> 

100.0 

100.0 

100.0 

ho  pyr.  <c+a> 

100.0 

100.0 

100.0 

7 

0.0023 

0.0023 

0.0023 

m 

0.02 

0.02 

0.02 

n 

0.01 

0.01 

0.01 

r 

0.1 

0.1 

0.1 

Table  3.5:  Calibration  of  Ti-3A1-2.5V  crystal  plast.  parameters. 


as  compared  with  experimental  results.  The  elastic  portion  of  these  curves  for  both 
RD  and  TD  specimens  matches  almost  perfectly.  Yield  strengths  also  look  good  for 
both  specimen  directions  with  the  RD  specimens  showing  a  lower  yield  point  consis¬ 
tent  with  experiments  and  the  TD  specimens  showing  a  slightly  higher  yield  point 
consistent  with  experiments.  Post-yield  slopes  leave  a  little  to  be  desired  when  it 
comes  to  matching  with  the  experiments  (at  least  for  1  of  the  2  directions).  For  MS- 
A  the  post-yield  slope  for  the  RD  specimen  look  almost  perfect;  however,  for  the  TD 
specimens  the  results  do  not  match  well.  Conversely,  for  MS-B  the  post-yield  slope 
for  the  TD  specimen  looks  good  while  the  RD  specimen  result  is  not  as  good.  It  is 
observed  that  after  about  1-3%  strain  that  the  post-yield  slopes  for  specimens  of  both 
orientations  converge  to  the  same  value  in  the  simulations,  while  in  the  experiments 
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this  is  not  the  case.  The  RD  and  TD  specimens  continue  to  show  a  marked  slope  dif¬ 
ference  past  1-3%  strain.  An  attempt  was  made  to  match  this  anisotropic  post-yield 
response  by  calibration  of  individual  plastic  slip  system  parameters.  Though  good 
agreement  could  be  achieved  by  modifying  ho  values  by  72  times  their  original  values, 
such  a  large  adjustment  of  previously  calibrated  parameters  did  not  seem  safe.  Fur¬ 
ther  study  is  recommended.  Particularly  needful  would  be  to  look  into  the  physical 
mechanisms  which  give  rise  to  such  a  dramatic  difference  in  post-yield  slopes  for  the 
2  specimen  orientations.  It  is  probable  that  the  current  models  are  not  able  (with¬ 
out  a  extremely  dramatic  adjustment  of  plastic  parameters)  to  capture  the  highly 
anisotropic  post-yield  slope  difference  past  1-3%  strains.  It  is  also  probable  that  error 
has  been  introduced  because  the  grain  size  validation  was  conducted  based  on  visual 
inspection  of  the  grain  size  distribution  based  on  number  fraction  instead  of  volume 
fraction.  A  further  cause  for  error  could  be  in  the  ODF  matching  (as  seen  in  the 
PFPD  plot  of  Figure  3.13)  and  sensitivities  to  perturbations  in  the  ODF  should  be 
studied  to  determine  if  indeed  the  ODF  matching  needs  to  be  improved. 

3.3.2  Model  Validation  Using  Load  Controlled  Experimental 
Data 

Model  validation  is  achieved  by  using  the  calibrated  model  to  predict  the  creep 
response  of  the  material  under  different  load  levels  and  specimen  orientations  for  MS - 
B.  Figure  3.22  shows  the  experimental  results  as  compared  to  the  simulation  results 
for  MS-B. 

Dwell  results  do  not  match  well  with  the  experimental  data  and  are  not  shown. 
The  current  model  has  no  parameter  or  equation  in  place  to  capture  the  mechanism 
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Figure  3.20:  Comparison  of  experimental  tensile  results  with  constant-strain  rate 
simulation  results  for  MS- A  (full  results  (a)  and  a  close-up  (b))  where  0°  strain  gage 
corresponds  to  i=j =1,  90°  strain  gage  corresponds  to  i=j =2,  and  45°  strain  gage 
corresponds  to  i=l,  j=2. 
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Figure  3.21:  Comparison  of  experimental  tensile  results  with  constant-strain  rate 
simulation  results  for  MS-B  (full  results  (a)  and  a  close-up  (b))  where  0°  strain  gage 
corresponds  to  i=j =1,  90°  strain  gage  corresponds  to  i=j =2,  and  45°  strain  gage 
corresponds  to  i=l,  j=2. 
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Figure  3.22:  Creep  validation  for  (a)  0°  strain  gages  and  (b)  90°  strain  gages. 


79 


that  gives  rise  to  such  a  dramatic  difference  in  plastic  strain  for  the  dwell  loading  ver¬ 
sus  the  creep  loading  (as  seen  in  Figure  3.8).  It  is  conjectured  that  the  reason  the  dwell 
experiments  have  such  a  high  plastic  strain  as  compared  to  the  creep  experiments  of 
the  same  maximum  load  is  due  to  dislocation  motion  during  the  unloading  and  reload¬ 
ing  portion  of  the  dwell  cycle.  During  unloading  and  reloading,  it  is  probable  that 
there  is  a  relaxation  effect  where  dislocations  piled  up  next  to  grain  boundaries  tend 
to  diffuse  into  the  center  of  grains  hireling  a  minimum  energy  configuration.  This 
corresponds  to  the  grains  being  able  to  see  more  plasticity  upon  the  next  dwell  hold. 
This  is  a  possible  reason  for  the  added  plastic  strain  for  the  dwell  loading  experiments 
as  compared  to  the  creep  loading  experiments.  The  current  model  does  not  include 
any  dependence  on  an  unloading/reloading  cycle  that  would  allow  for  this  relaxation 
and  the  associated  softening.  This  could  be  included  in  future  work  which  could  take 
advantage  of  the  large  amount  of  dwell  experimental  data  presented  here. 
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Chapter  4:  Conclusions  and  Future  Work 


This  thesis  demonstrates  that  a  complete  procedure  is  in  place  to  study  a  large 
number  of  material  microstructures  by  computational  means.  A  relatively  efficient 
data  collection  and  processing  procedure  and  a  robust  synthetic  microstructure  gen¬ 
erating  program  combined  with  a  meshing  software  and  crystal  plasticity  based  finite 
element  analysis  program  were  integrated  to  yield  a  useful  computational  tool.  This 
procedure  is  able  to  accurately  predict  the  macroscopic  constant  strain  rate  response 
as  well  as  the  primary  and  secondary  macroscopic  creep  response  of  important  high- 
performance  titanium  alloys.  Because  the  model  calculates  the  elastic-plastic  response 
of  individual  grains  along  specific  slip  systems,  insight  is  gained  into  the  ”  inner  work¬ 
ings”  of  the  material  response  on  a  grain  level.  In  addition,  large  sets  of  model 
microstructures  with  vastly  different  characteristics  were  created  and  used  to  develop 
functional  forms  that  relate  microstructural  parameters  to  responses  such  as  yield 
strength,  creep  parameters,  and  Ramberg-Osgood  parameters.  The  proposed  model 
for  yield  strength  as  a  function  of  a  lath  thickness  is  validated  against  a  previous 
published  regression  model  based  off  75  separate  experiments. 
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In  addition,  a  CPFEM  model  update  has  been  proposed  and  implemented  that 
provides  strides  towards  modeling  plasticity  more  realistically  in  grains  with  elon¬ 
gated  shapes.  Because  of  the  large  amount  of  experimental  data  available  for  speci¬ 
mens  of  different  orientations,  the  current  models  were  able  to  be  calibrated  as  well  as 
validated  by  predicting  the  creep  response  under  different  loads  and  specimen  orien¬ 
tations.  The  current  model  is  not  able  to  capture  the  dwell  fatigue  response;  however, 
the  experimental  data  reported  here  could  be  used  in  the  future  to  develop  a  model 
that  captures  this  pronounced  dwell  effect. 

Future  work  could  also  consider  time  to  failure  under  fatigue  loading  as  a  response 
parameter.  The  ” weakest  link,”  local  nature  of  fatigue  failure  would  require  the 
introduction  of  probabilistic  methods  and  the  large  times  to  failure  of  many  materials 
requires  multi-scaling  techniques  but  nonetheless  the  groundwork  is  laid  for  developing 
these  important  functional  dependencies. 

In  addition,  the  functional  dependencies  developed  in  this  work  are  a  step  towards 
helping  to  remove  the  industry’s  heavy  reliance  on  expensive  experiments  as  well  as 
their  need  to  perform  expensive  crystal  plasticity  simulations.  The  sensitivity  analyses 
performed  here  also  give  feedback  to  those  responsible  for  advancing  the  accuracy  and 
efficiency  of  materials  characterization  as  to  the  microstructural  features  that  are  the 
most  critical  to  particular  mechanical  responses.  This  information  can  be  valuable  to 
the  end  that  time  and  energy  can  be  focused  on  advancing  the  efficiency  and  accuracy 
of  characterization  of  the  features  most  crucial  to  particular  response  variables. 
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